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ABSTRACT 

The  more-complete-expansion  cycle  can  be  applied 
to  a  stationary  irrigation  engine  commonly  operated  at 
constant  partial-load  burning  liquid  fuels.   The 
analysis  involved  the  development  of  a  digital 
computer  model  to  simulate  a  real  engine. 

The  Chevrolet  5.74  liter  (350  cu.  in.)  engine  was 
simulated  to  determine  the  performance  changes  and 
possible  advantages  of  applying  the  more-complete- 
expansion  cycle  features.   The  fuel  used  in  this  model 
was  methane  gas.   The  generated  performance  data  is 
presented  in  figures.   The  modified  engine  specifications 
are  results  of  these  figures  and  are  used  to  determine 
the  cost  effectiveness  of  the  concept. 

The  performance  results  in  conjunction  with  a 
cost  analysis  suggest   that  a  fuel  saving  greater  than 
11%  can  be  achieved  by  utilizing  the  more-complete- 
expansion  cycle.   If  the  engine  is  operated  frequently, 
the  amount  of  money  saved  in  fuel  is  significant  and 
could  pay  off  the  cost  of  the  engine  modification  in 
less  than  a  year.   Therefore,  this  idea  is  recommended 
for  further  development  and  testing. 
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1.0  INTRODUCTION 

The  problem  exists  that  some  agricultural  locations 
along  the  "Ogallala  Aquifer"  use  ordinary  automobile  or 
tractor  spark  ignition  engines  to  run  irrigation  pumps. 
These  engines  generally  burn  natural  gas  (about  95% 
methane)  and  operate  at  part-throttle,  using  only  a 
fraction  of  the  useful  displacement  volume  while 
consuming  large  amounts  of  fuel.   It  is  proposed  that 
the  engine  can  be  modified  in  such  a  way  to  incorporate 
a  more-complete-expansion  process.   This  process  involves 
the  use  of  an  expansion  ratio  greater  than  the  effective 
compression  ratio.   The  addition  of  the  more-complete- 
expansion  process  tends  to  use  the  entire  displacement 
volume  more  efficiently,  thus  decreasing  the  fuel 
consumption . 

The  purpose  of  this  research  is  to  determine  the 
merit  of  implementing  such  a  process  in  a  stationary 
internal  combustion  engine.   The  engine  will  be  a 
Chevrolet  5.74  liter  V-8  (350  cu.  in.).   A  computer 
program  will  be  used  to  simulate  the  part-throttle 
reference  case  and  the  modified  cases  to  determine  the 
advantages,  if  any. 
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2.0  LITERATURE  REVIEW 

The  more-complete-expansion  cycle  known  as  the 
Atkinson  Cycle  was  named  after  Atkinson  who  first 
proposed  the  idea  in  1885.   This  cycle  has  not  been 
researched  or  used  significantly  in  the  past.   Only 
three  references  found  mention  any  work  or  use  of  such 
a  cycle.   The  reason  for  this  lack  of  interest  is  due 
to  the  fact  that  the  more-complete-expansion  cycle  is 
most  effective  when  the  engine  is  operated  under  a 
constant  heavy  load.   Another  disadvantage  of 
utilizing  this  cycle  is  the  size  of  the  engine  required 
For  this  cycle  to  be  successful,  a  large  displacement 
volume  is  needed  to  compensate  for  the  lower  mean 
effective  pressure  which  results. 

The  first  reference  is  a  book  entitled,  "The 
Internal  Combustion  Engine  in  Theory  and  Practice", 
Vol.  II  (p. 402),  by  Charles  Fayette  Taylor.   This 
reference  includes  an  ideal  theoretical  example  of 
designing  a  large  gas  engine  to  operate  a  natural 
gas  compressor  unit.   The  purpose  of  using  the  more- 
complete-expansion  cycle  is  three-fold.   First,  the 
compressor  station  requires  great  reliability,  second, 
long  life  and  third,  good  fuel  economy.   All  these 
characteristics  can  be  achieved  by  utilizing  the  more- 
complete-expansion  cycle. 
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The  second  reference  is  a  technical  paper 
entitled,  "Improved  Engine  Efficiency  with  Emphasis 
on  Expansion  Ratios",  by  Rodney  Herbert  and  brian 
Lanoway.   The  paper  discusses  ways  in  which  the 
expansion  ratio  effects  the  efficiency  of  an  engine. 
It  mentions  the  use  of  a  modified  Atkinson  Cycle  where 
the  intake  valve  is  closed  late.   This  allows  some  of 
the  freshly  inducted  charge  to  be  forced  back  into  the 
intake  manifold  which  effectively  reduces  the 
compression  ratio. 

The  third  is  a  book  entitled,  "Combustion  Engine 
Processes",  by  Lester  C.  Lichty.   The  information 
included  in  this  reference  deals  with  the  ideal 
advantages  of  applying  the  more-complete-expansion 
cycle  . 

Late  valve  timing  has  been  examined  by  Miller 
with  the  Cadillac  and  Nordberg  engines  and  lately  by 
SAAB  and  Edward  LaForce.  (Herbert,  Lanoway,  1979,  pl5) 
LaForce  claimed  that  his  engine  operating  on  the 
modified  Atkinson  Cycle  demonstrated  a  25%  increase  in 
fuel  economy,  however  the  American  EPA  tested  the 
engine  and  found  it  to  give  better  fuel  economy  but 
only  at  the  expense  of  reduced  power  output  and 
increased  emissions.   The  paper  also  suggests  that 
SAAB  found  a  slight  increase  in  fuel  economy  but  feels 
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the  cost  of  production  is  unwarranted. 

Not  having  found  any  useful  references  through 
published  literature  or  the  library  computerized 
search  in  regards  to  the  more-complete-expansion 
cycle,  this  study  will  develop  a  computer  simulation 
of  a  more-complete-expansion  cycle. 
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3.0  PROGRAM  DEVELOPMENT 
3.1  Description: 

The  computer  program  listed  in  Appendix  III 
used  in  this  research  was  developed  largely  from 
the  discussion  in  Ashley  S.  Campbell's  book 
entitled,  "Thermodynamic  Analysis  of  Combustion 
Engines".   The  program  simulates  the  processes  of 
an  actual  engine  and  generates  performance 
information.   Some  of  the  information  produced  by 
the  program  includes  the  spark  advance,  mean  piston 
speed,  mean  effective  pressure,  brake  power  output, 
brake  specific  fuel  consumption,  volumetric 
efficiency,  pressure  before  the  compression  stroke, 
peak  cylinder  pressure  and  brake  thermal  efficiency. 

The  equations  used  to  calculate  the  performance 
information  are  shown  below. 

MEP  -   WNET 


EFF  = 


VDISP 

-100*WNET*NM0 
NM-URP 


POWER  =  WNET»RPM*f^l*fc  CyCle] 

LoOsecJ  Lz  rev   J 


SFf  =  — - — 

POWER 
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where  , 

WNET=  net  work  per  cycle,  kJ 

VDISP=  displacement  volume,  m3 

NMO=  stoichiometric  moles,  kmol 

NM=  number  of  kmol  of  fuel  mixture,  kmol 

URP=  constant  volume  heat  of  reaction,  kJ/kmol 

RPM=  engine  speed,  rev/min 

FF=  fuel  flow,  g/hr 

POWER=  engine  power  per  cycle,  kW 

MEP=  mean  effective  pressure,  kPa 

EFF=  brake  thermal  efficiency,  % 

SFC=  brake  specific  fuel  consumption,  g/kW-hr 

The  model  takes  into  account  progressive 
combustion,  realistic  valving  and  heat  transfer 
during  the  compression,  combustion  and  expansion 
processes.   Friction  and  gas  properties  are 
determined  from  correlated  functions.   One  of  the 
largest  uncertainties  in  the  program  is  the  heat 
transfer  to  or  from  the  cylinder  walls. 
Specifically,  there  exists  no  precise  approximation 
to  the  heat  transfer  coefficient  inside  the  cylinder 
during  the  cycle.   This  is  an  area  of  present 
research.   The  program  does  include  a  correlation 
for  the  heat  transfer  coefficient.   This  correlation 
appears  to  produce  satisfactory  coefficients  for 
use  in  this  model. 

As  noted  before,  the  computer  model  was 

developed  primarily  from  Campbell's  book.   However, 

there  are  three  important  features  included  in  the 

program  which  are  not  discussed  completely  in  this 

text.   The  features  are  heat  transfer,  mechanical 
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friction  and  mass  flow  coefficients.   Each  of  these 
areas  were  added  through  research  of  published 
literature  . 

The  instantaneous  heat  transfer  coefficient 
correlation  was  found  in  a  technical  paper  entitled, 
"Convective  Heat  Transfer  in  Reciprocating  Engines", 
by  Bardon  and  Rao.   The  paper  uses  flow 
visualization  studies  and  probe  explorations  to 
determine  the  complex  nature  of  the  turbulent  flow 
inside  a  cylinder  of  an  internal  combustion  engine. 
These  results  were  combined  with  convection  heat 
transfer  theory  to  produce  the  following  correlation. 


H=.0006* 


S*CP*P 
T*R 


*[¥] 


where  , 

H=  inst  .  heat  transfer  coefficient,  kW/m2-K 

S=  stroke ,  m 

CP=  constant  pressure  specific  heat,  kJ/kmol-K 

P=  inst.  pressure,  kPa 

VIC=  cylinder  volume  at  intake  closing,  m3 

V=  inst.  volume,  m3 

T=  inst.  temperature,  K 

R=  universal  gas  constant,  kJ/kmol-K 

This  correlation  produced  instantaneous  heat 
transfer  coefficients  of  200-3000  W/m2-K.   The 
occurrence  of  these  values  during  the  cycle  have  been 
examined  and  proven  satisfactory  for  use  in  the  model 

Another  area  of  improvement  in  the  model  is  the 
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discharge  coefficients  during  the  valving  events. 
The  correlations  used  to  determine  the  flow 
coefficients  was  taken  from  Taylor's  text  (Vol.  I, 
p. 173).   It  is  listed  below. 


•Is-  -3]    +Hi"  -3]  +  -7 


-20 


where  , 

CD=drag  coefficient 

L  =  valve  lif t ,  m 

D=  valve  diameter,  m 

The  drag  coefficient  is  a  function  of  the  lift/valve 
diameter  ratiowhich  is  found  instantaneously  during  the 
valving  subroutines. 

The  final  feature  added  is  a  simple  correlation 
for  engine  friction  in  a  standard  4-stroke  spark 
ignition  engine.   This  correlation  shown  below  was 
also  found  in  Taylor's  text  (Vol.  II,  p. 382). 

FMEP=  ^p-*   [.076*  [MPS]  2  +  .92*  [MPS]  +  7.23J 

where  , 

FMEP=frictional  mean  effective  pressure,  kPa 
MPS  =  mean  piston  speed,  m/s 

Figures  la  &  lb  show  the  effect  of  adding  the 

three  non-ideal  conditions  to  the  model  derived  from 

Campbell's  book.   The  simulations  were  conducted  on  a 

Chevy  5.74  liter  engine  burning  gasoline.   Curve  one 

is  the  model  developed  from  the  text.   The  second 

curve  is  the  same  model  with  the  addition  of  flow 


19 


210 

200 
190 

isa 

170 

160 
150 

140 

130 

120 

^  no 

6  100 

u 

0) 

§    90 

Oh 

g    804 


00 

c 
w 


70 

60- 

50- 

40 

30- 

20 

10 


1  -  Text 

2  -  Flow  Coefficient 

3  -  Heat  Transfer 

4  -  Friction 


2000 


3000 
Engine  Speed  (RPM) 


4000 


Figure  la-  Engine  power  versus  engine  speed  for  a  Chevy  5.74  liter 
engine  with  various  non-ideal  processes. 
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Figure  lb-  Brake  specific  fuel  consumption  versus  engine 
speed  for  a  5.74  liter  engine  with  various 
non-ideal  processes. 
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coefficients  during  valving.   The  third  curve  shows  the 
characteristics  of  the  model  including  flow  coefficients 
and  heat  transfer.   The  final  curve,  number  4,  which 
best  approximates  the  actual  engine,  is  the  model  used 
in  this  work  including  flow  coefficients,  heat  transfer 
and  friction . 

Other  irreversible  processes  such  as  inertial 
effects,  chemical  equilibrium,  heat  transfer  during 
valving,  pumping  losses  and  parasitic  losses  are  not 
included  in  this  model.   These  losses  combine  for  at 
least  a  10%  uncertainty  within  the  program. 

The  input  parameters  for  the  model  are  the 
compression  ratio,  type  of  fuel,  air/fuel  ratio,  bore 
(cm),  stroke  (cm),  connecting  rod  length  (cm),  valve 
lift  (cm),  valve  dimensions  (cm),  valve  timing,  exhaust 
pressure  (kPa),  engine  speed  (RPM),  intake  manifold 
pressures  (kPa)  and  temperatures  (K).   All  of  these 
parameters  are  important  characteristics  necessary  for 
an  approximate  simulation  of  an  actual  engine. 
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3.2  Validity: 

After  developing  the  model  it  was  necessary  to 
determine  the  validity  of  its  predictions.   Since  all 
engines  can  be  stocked  or  modified  in  several  different 
configurations,  it  was  difficult  to  receive  performance 
information  from  a  professional  testing  agency  for  the 
particular  engine  in  question.   The  predicted 
performance  information  for  a  Chevy  5.74  liter  V-8  engine 
burning  gasoline  was  compared  with  some  good  experimental 
data  which  Dr.  Don  Fenton,  professor  at  Kansas  State 
University,  had  taken  while  employed  at  New  Mexico  State 
University.   The  results  of  the  comparison  can  be  seen 
in  Figures  2  &  3.   From  Figure  2  note  that  the  computer 
simulation  predicts  the  experimental  data  to  within 
about  10  percent.   Agreement  is  better  at  high  engine 
speeds  --  above  3800  RPM.   Note  also  that  the  uncertainty 
is  consistently  positive  (experimental  data  less  than 
predicted  power)  and  so  the  model  is  expected  to  yield 
accurate  trends  associated  with  comparisons.   Figure  3 
suggests  an  uncertainty  of  less  than  8%  between  model 
prediction  and  experimental  data  for  the  BSFC.   Therefore, 
it  is  assumed  that  the  model  exhibits  enough  validity  to 
continue  with  the  modeling  of  a  more-complete-expansion 
cycle . 
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4.0  DISCUSSION  OF  REFERENCE  CASE 
The  reference  engine  will  be  a  Chevy  5.74  liter 
V-8  operated  at  part-throttle,  approximately  178  mm  Hg 
vacuum,  and  standard  atmospheric  conditions.   The 
reason  for  selecting  this  type  of  internal  combustion 
engine  is  due  to  its  popularity  and  accessibility  of 
modified  parts.   The  input  parameters  for  the  reference 
case  are  listed  in  Appendix  I.   These  values  are  true 
measurements  for  a  stock  Chevy  5.74  liter  V-8  engine. 

Simulations  were  conducted  on  the  reference  engine 
to  determine  the  correct  air/fuel  ratio  for  maximum 
performance.   The  stoichiometric  air/fuel  ratio  for 
methane  gas  is  2.   This  means  that  for  every  1  mole  of 
fuel  2  moles  of  air  are  required  for  complete  combustion 
Complete  combustion  refers  to  the  condition  where  the 
products  contain  neither  carbon  monoxide  nor  excess  air. 
This  condition  also  results  in  the  greatest  adiabatic 
flame  temperature  which  is  a  very  important  factor  in  an 
internal  combustion  engine. 

Experimental  tests  show  that  maximum  power  output 
and  fuel  efficiency  are  achieved  at  the  stoichiometric 
condition  or  when  slightly  lean.  (Ferguson,  1986). 
Figure  4  shows  the  results  from  the  simulation.   It 
shows  that  the  model  predicts  maximum  performance 
occurring  at  a  lean  air/fuel  ratio  of  2.2.   This  value 
will  be  used  in  the  remaining  simulations. 
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Table  1  is  a  list  of  performance  data  generated 
by  the  computer  for  the  reference  case  varying  the  RPM 
The  information  shown  will  be  useful  for  comparison 
purposes  as  the  research  progresses. 
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5.0  MORE-COMPLETE-EXPANSION  CYCLE 
The  more-complete-expansion  cycle  is  analogous  to 
a  standard  Otto  Cycle  except  for  the  increased 
expansion  ratio  and  decreased  compression  ratio.   Both 
cycles  involve  compression,  combustion,  expansion, 
exhaust,  and  intake  of  a  certain  fuel  which  is  burnt  to 
produce  useful  power  output.   In  the  more-complete- 
expansion  cycle,  the  ignition  pressure  is  controlled  by 
early  or  late  closing  of  the  intake  valve.   Typically, 
the  cycle  in  which  the  intake  valve  is  closed  early  is 
known  as  an  Atkinson  Cycle  and  the  other  is  called  a 
Modified  Atkinson  Cycle.   Whichever  cycle  is  chosen  the 
results  are  usually  very  similar. 

The  expansion  process  in  an  ideal  more-complete- 
expansion  cycle  is  lengthened  by  late  opening  of  the 
exhaust  port,  thus  allowing  expansion  of  the  gas  to 
continue  until  it  reaches  atmospheric  pressure.   Upon 
reaching  atmospheric  pressure,  the  piston  is  usually  at 
bottom  dead  center  at  which  time  the  exhaust  valve  opens 
and  the  piston  pushes  the  gases  out  of  the  cylinder. 
However,  this  type  of  cycle  is  not  viable  for  engines 
that  operate  at  partial-load  because  in  light  load 
situations  the  mean  cylinder  pressure  at  the  end  of 
expansion  tends  to  be  near  or  below  the  pressure 
required  to  overcome  friction.   Under  such  circumstances 
the  application  of  a  more-complete-expansion  process 
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leads  to  a  net  loss  rather  than  a  gain  in  efficiency. 

The  reason  for  this  arrangement  is  that  maximum 
cylinder  pressure  is  usually  limited  by  detonation  or 
knock.   In  the  more-complete-expansion  cycle, 
controlling  the  ignition  pressure  results  in  a  more 
controllable  maximum  cylinder  pressure.   In  the  past, 
spark  ignition  engines  used  increased  compression 
ratios  to  increase  cycle  efficiency.   This  was 
severely  restricted  because  of  detonation  and 
requirements  for  higher  octane  fuels.   Therefore,  the 
only  other  feasible  way  to  significantly  increase  cycle 
efficiency  is  to  increase  the  expansion  ratio. 

In  the  simulations  to  follow,  this  investigation 
will  attempt  to  determine  the  effects  of  adding  such  a 
process  to  an  irrigation  engine  that  operated  at 
part-throttle.   The  results  will  be  compared  to  the 
reference  case  to  determine  the  advantages  of  the 
modification . 
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5.1  Matched  ignition  pressure  with  late  intake  closing: 
The  technique  used  in  the  simulations  was  first 
to  open  the  throttle  plate  wide  open  allowing  the 
intake  manifold  pressure  to  be  atmospheric.   Secondly, 
while  increasing  the  compression  ratio  the  intake 
valve  closing  was  delayed  until  the  pressure  prior 
to  ignition  was  the  same  as  the  ignition  pressure 
in  the  reference  engine.   The  reason  for  matching 
the  ignition  pressures  is  so  all  the  performance 
information  generated  by  the  model  will  have  one 
point  in  common.   Comparing  the  two  results  provide 
a  useful  way  to  determine  how  the  more-complete- 
expansion  process  effects  the  remaining  cycle. 

The  increase  in  compression  ratio  not  only 
improves  efficiency  but  also  helps  offset  the 
decrease  in  peak  cylinder  pressure  caused  by  the 
late  closing  of  the  intake  valve. 

The  results  of  the  first  simulation  varying  the 
compression  ratio  from  8.5:1  to  13.5:1  while  matching 
ignition  pressures  can  be  seen  in  Figures  5a-5g. 
Figure  5a  illustrates  how  the  compression  ratio 
affects  efficiency.   It  shows  that  as  the  compression 
ratio  increases  the  efficiency  increases,  but  only  to 
a  point  where  it  levels  off,  namely  at  a  compression 
ratio  of  about  12.5:1.   This  could  have  been 
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Figure  5a  -  Efficiency  versus  compression  ratio  for 
a  Chevy  5.74  liter  with  matched  ignition 
pressures,  (late  intake  valve  closing) 
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gure  5b  -  Engine  power  versus  compression  ratio 
for  a  Chevy  5.74  liter  with  matched 
ignition  pressures,  (late  intake  valve 
closing  ) 
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Figure  5c  -  Brake  specific  fuel  consumption  versus  compression 
ratio  for  a  Chevy  5.74  liter  with  matched  ignition 
pressures,  (late  intake  valve  closing) 
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predicted.  Similar  results  occur  for  an  air-standard 
cycle,  except  for  the  rapidly  leveling  off  of  the 
curve . 

Figure  5b  shows  the  effect  of  the  compression 
ratio  on  the  power  output.   This  curve  continues  to 
decrease  as  the  compression  ratio  increases.   The 
decrease  in  power  output  is  due  to  the  fact  that  the 
intake  valve  closing  is  continually  retarded  in 
order  to  match  the  ignition  pressures.   This 
retarding  causes  more  of  the  freshly  ingested  charge 
to  be  expelled  back  into  the  intake  manifold,  thus 
decreasing  the  amount  of  fuel  available  for 
conversion  to  work.   The  exact  effect  of  pushing 
some  of  the  heated  charge  into  the  intake  manifold 
is  not  known,  however,  increasing  the  intake 
manifold  temperature  by  8%  tends  to  decrease  the 
thermal  efficiency  by  4%  and  the  power  output  by 
12.3%.   This  situation  is  beyond  the  scope  of  this 
research  because  of  the  complexity  involved  in 
determining  the  correct  intake  manifold  temperature. 
Should  this  configuration  be  chosen  for 
experimentation,  this  effect  must  be  examined  since 
it  tends  to  deteriorate  engine  performance. 

Figure  5c  is  a  plot  of  brake  specific  fuel 
consumption  versus  compression  ratio.   This  curve  is 
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inversely  similar  to  the  efficiency  curve  since 
they  are  essentially  dependent  upon  one  another. 
The  minimum  BSFC  point  occurs  at  a  compression  ratio 
of  12.5:1. 

As  can  be  seen  in  Figure  5d ,  the  fuel  flow 
behaves  like  the  power  curve.   Since  the  compression 
ratio  is  increased  it  is  necessary  to  retard  the 
intake  valve  closing  so  the  pressure  at  ignition  is 
matched.   In  the  process  of  closing  the  intake  valve 
later,  more  of  the  fresh  charge  is  pushed  out  of  the 
cylinder  causing  a  reduction  in  fuel  flow. 

An  unexpected  situation  arises  when  the  peak 
pressure  is  plotted  versus  compression  ratio.   One 
would  expect  the  peak  pressure  to  be  lower  since  a 
smaller  amount  of  fuel  is  being  burned  as  the 
compression  ratio  increases.   This  is  not  the  case. 
(See  Figure  5e)   In  fact,  the  peak  pressure 
increases  linearly  as  the  compression  ratio 
increases . 

The  underlying  causes  of  this  behavior  are  not 
completely  understood.   The  analysis  of  this 
condition  is  complex  and  only  a  speculation  can  be 
made  to  determine  the  reasons  for  the  linear 
behavior.   Since  the  only  two  parameters  which 
changed  from  cycle  to  cycle  were  the  compression 
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ratio  and  valve  timing,  it  is  understood  that  the 
peak  cylinder  pressure  is  a  function  of  both 
parameters.   The  shape  of  the  actual  function  may 
or  may  not  be  linear.   The  reason  for  this 
uncertainty  is  because  of  the  simple  model  within 
the  program  used  for  the  simulation  of  the  combustion 
process.   If  a  more  accurate  function  is  sought,  a 
more  detailed  and  complex  model  will  be  necessary. 

It  has  been  determined  that  the  two  factors 
involved  in  this  behavior  are  the  change  in  valve 
timing  and  an  increased  compression  ratio.   If  each 
parameter  is  examined  individually  a  better 
understanding  of  the  speculation  will  be  achieved. 

The  result  of  retarding  the  valve  timing  while 
holding  the  compression  ratio  constant  is  a  reduced 
number  of  moles  of  fuel  mixture  ingested  into  the 
cylinder  to  be  burned.   The  peak  cylinder  pressure 
and  the  mean  effective  pressure  are  lowered  due  to 
a  smaller  amount  of  fuel  available  for  work. 

The  result  of  increasing  the  compression  ratio 
while  holding  the  number  of  moles  of  fuel  constant 
is  an  increased  peak  cylinder  pressure.   This  is 
caused  by  a  change  in  the  cylinder  pressure  during 
combustion  which  is  broken  into  two  parts  - 
compression  and  combustion.   Since  the  amount  of 
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fuel  is  constant,  the  pressure  change  associated  with 
combustion  regains  the  same.   Therefore,  the 
increased  compression  ratio  causes  the  compressive 
pressure  change  to  be  larger.   This  results  in  a 
greater  peak  cylinder  pressure. 

When  the  two  conditions  occur  simultaneously, 
each  of  the  individual  behaviors  combine  and  offset 
one  another  to  produce  the  linear  function. 
Remember,  each  of  these  functions  are  unknown  and 
it  is  only  speculated  that  they  combine  to  produce 
the  straight  line.   As  stated  earlier,  a  detailed 
examination  of  each  individual  function  is  needed 
in  order  for  a  complete  understanding  of  this 
behavior  . 

The  following  two  figures  (Figures  5f  ■  5g)  are 
curves  demonstrating  how  the  more-complete-expansion 
cycle  effects  the  mean  effective  pressure  and  the 
volumetric  efficiency.   The  mean  effective  pressure 
remains  somewhat  constant  while  the  volumetric 
efficiency  deteriorates  due  to  the  retarding  of  the 
intake  valve  . 

Another  concern  regarding  application  of  the 
more-complete-expansion  cycle  is  which  exhaust 
valve  opening  position  will  maximize  the  cycle 
efficiency.   From  the  previous  analysis,  it  appears 
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that  increasing  the  compression  ratio  above  12.5:1 
is  unwarranted  since  cycle  efficiency  tends  to 
level  out.   It  is  unwarranted  because  of  the 
additional  peak  stresses  placed  on  the  cylinder  to 
produce  only  a  small  change  in  efficiency. 

For  the  purpose  of  determining  the  optimum 
exhaust  valve  timing  a  compression  ratio  of  12.5:1 
was  chosen.   The  results  from  the  model  are  shown  in 
Figures  6a-6c .   All  of  these  figures  show  that  the 
optimum  exhaust  valve  opening  is  20  degrees  BBDC. 
The  reason  for  opening  the  exhaust  valve  before 
bottom  dead  center  is  because  the  remaining  pressure 
in  the  cylinder  is  needed  to  overcome  the  drag 
forces  on  the  exhaust  gases  caused  by  the  valve 
dynamics . 

On  the  basis  of  the  previous  simulations,  the 
model  predicts  that  the  greatest  advantage  in 
performance  can  be  achieved  by  modifying  the 
reference  engine  in  the  following  ways.   Increase 
the  compression  ratio  to  12.5:1,  intake  valve  closinj 
as  shown  in  Table  2  and  the  exhaust  valve  opening  at 
20  degrees  BBDC.   The  model  predictions  for  this 
configuration  are  shown  in  Table  2  and  the  P-V 
diagrams  of  this  cycle  and  the  reference  cycle  are 
shown  in  Figure  7. 
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Figure  7 
P-V  diagrams  of  Base  Cycle  and 
Modified  Cycle  having  matched 
ignition  pressure  and  late 
intake  valve  closing,  2000  RPM 
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5.2  Matched  ignition  pressure  early  intake  closing: 

The  technique  used  in  the  following  simulation 
is  similar  to  those  used  in  the  late  closing  case. 
In  the  early  intake  closing  case  the  intake  valve 
is  closed  before  the  intake  process  is  complete. 
The  gases  in  the  cylinder  are  then  pulled  below 
atmospheric  pressure  until  bottom  dead  center.   After 
bottom  dead  center,  the  piston  begins  to  compress  the 
gases  until  the  appropriate  ignition  pressure  is 
reached.   The  advantage  in  using  this  type  of 
arrangement  is  that  none  of  the  heated  gases  are 
forced  back  into  the  intake  manifold  causing  a 
temperature  increase.   The  purpose  of  closing  the 
intake  valve  early  and  then  pulling  the  gases  in  the 
cylinder  toward  a  vacuum  is  to  effectively  reduce 
the  compression  ratio.   The  cycle  analysis  of  this 
situation  is  much  the  same  as  before. 

Once  again  the  ignition  pressures  are  matched, 
compression  ratio  varied  and  the  intake  valve  is 
closed  early.   In  Figures  8a-8g  the  effects  of  the 
more-complete-expansion  cycle  and  varying  the 
compression  ratio  can  be  seen.   Each  of  the  curves 
behave  similarly  to  the  case  in  which  the  intake 
valve  closed  late.   Therefore,  the  best  compression 
ratio  to  select  without  stress  consideration  is 
12.5:1. 
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Figure  8b  -  Engine  power  versus  compression  ratio  for  a 

Chevy  5.74  liter  with  matched  ignition  pressures, 
(early  intake  valve  closing) 
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Figure  8c  -  Brake  specific  fuel  consumption  versus  compression 
ratio  for  a  Chevy  5.74  liter  with  matched  ignition 
pressures,  (early  intake  valve  closing) 
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Figure  8d  -  Fuel  flow  versus  compression  ratio  for  a  Chevy  5.74 
liter  with  matched  ignition  pressures,  (early 
intake  valve  closing) 
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Figure  8e  -  Peak  cylinder  pressure  versus  compression  ratio 
for  a  Chevy  5.74  liter  with  matched  ignition 
pressures,  (early  intake  valve  closing) 
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Figure  8f  -  Volumetric  efficiency  versus  compression 

ratio  for  a  Chevy  5.74  liter  with  matched 
ignition  pressures,  (early  intake  valve  closing) 
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Figure  8g  -  Mean  effective  pressure   versus  compression  ratio 
for  Chevy  5.74   liter      with  matched   ignition 
pressures,    (early   intake   valve  closing) 
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Figures  9a-9c  are  plots  to  determine  the 
optimum  exhaust  valve  opening.   As  before,  the 
optimum  time  is  20  degrees  BBDC.   Table  3  shows  the 
results  of  the  proposed  arrangement  with  early 
intake  closing.   Figure  10  is  a  plot  of  the  P-V 
diagram  of  this  cycle  compared  to  the  reference 
cycle . 
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Figure  9b  -  Engine  power  versus  exhaust  opening  for  a 
Chevy  5.74  liter  with  CR=12.5  and  matched 
ignition  pressures,  (early  intake  valve  closing) 
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Figure  10 

P-V  diagrams  of   Base  Cycle  and 

Modified   Cycle   having  matched 

ignition   pressure  and   early 

intake   valve  closing,    2000  RPM 
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5.3  Matched  power  output,  late  intake  valve  closing: 

The  purpose  of  this  segment  in  the  analysis  is 
to  determine  the  advantage  of  implementing  the  more- 
complete-expansion  idea  to  an  actual  irrigation 
pump  arrangement.   Irrigation  pump  gear  heads 
usually  require  a  specified  amount  of  power  at  a 
certain  speed.   It  is  essential  that  the  modified 
engine  produce  the  required  pumping  conditions. 
These  conditions  are  met  by  matching  the  predicted 
power  output  from  the  modified  model  to  those  of 
the  reference  model  which  was  operated  at  178  mm  Hg . 
This  also  requires  that  the  engine  speed  (RPM)  be 
equivalent  for  an  accurate  prediction. 

The  procedure  used  for  matching  the  power 
output  is  much  the  same  as  the  two  previous  cases. 
The  match  was  achieved  by  varying  the  intake  valve 
closing  while  maintaining  the  exhaust  valve  opening 
at  its  optimum  time  of  20  degrees  BBDC. 

Simulations  were  produced  to  examine  the  effects 
of  changing  the  compression  ratio  and  intake  valving 
on  the  performance  information.   In  this  case  the 
intake  valve  was  closed  late  allowing  a  large 
portion  of  the  freshly  inducted  charge  to  be  pushed 
out  of  the  cylinder  resulting  in  a  reduced  power 
output.   The  results  of  the  simulations  can  be  seen 
in  Figures  lla-lle. 
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Figure  11a  -  Efficiency  versus  compression  ratio  for 
a  Chevy  5.74  liter  with  matched  power 
output,  (late  intake  valve  closing) 
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compression  ratio  for  a  Chevy  5.74  liter 
with  matched  power  output,  (late  intake 
valve  closing) 
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for  a  Chevy  5.74  liter   with  matched  Power  output 
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Figure  lie-  Volumetric  efficiency  versus  compression  ratio 

for  a  Chevy  5.74  liter   with -matched  power  output 
(late  intake  valve  closing) 
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From  these  figures  it  is  observed  that  a 
significant  increase  in  engine  performance  cannot  be 
attained  by  increasing  the  compression  ratio  beyond 
12.5:1.   This  condition  was  chosen  as  the  maximum 
performance  arrangement  for  matching  power  output. 
Table  A  shows  the  performance  data  for  a  simulated 
engine  with  matched  power  output,  late  intake  valve 
closing  and  12.5:1  compression  ratio.   The 
information  produced  will  be  useful  in  the  comparison 
between  the  modified  engine  and  the  reference  engine. 
Also,  a  P-V  diagram  of  the  modified  cycle  and 
reference  cycle  is  included  in  Figure  12  to 
demonstrate  the  changes  in  the  two  cycles. 
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Figure  12 
P-V  diagrams  of  Base  Cycle  and 
Modified  Cycle  having  matched 
power  output  (kW)  and  late 
intake  valve  closing,  2000  RPM 
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5.4  Matched  power  output,  early  intake  closing: 

This  section  involves  determining  the  effects 
on  engine  performance  when  closing  the  intake  valve 
early  while  holding  the  power  output  the  same  as  the 
reference  case.   The  method  used  in  this  case  is  the 
same  as  before  and  the  plotted  results  are  shown  in 
Figures  13a-13e.   These  figures  show  that  the  best 
compression  ratio  for  maximum  performance  under 
matched  power  output  constraints  is  12.5:1.   This 
value  is  used  to  produce  Table  5  which  lists  the 
performance  data  resulting  from  the  change.   A  P-V 
diagram  is  included  in  Figure  14  to  show  the 
modifications  between  the  reference  and  modified 
cases  . 
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Figure  13a  -  Efficiency  versus  compression  ratio  for 
a  Chevy  5.74  liter  with  matched  power 
output,  (early  intake  valve  closing) 
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Figure  13b  -  Brake  specific  fuel  consumption  versus 

compression  ratio  for  a  Chevy  5.74  liter 
with  matched  power  output,  (early  intake 
valve  closing) 
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Figure  13c  -  Fuel  flow  versus  compression  ratio  for 
Chevy  5.74  liter  with  matched  power 
output,  (early  intake  valve  closing) 
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Figure  13d-  Peak  cylinder  pressure  versus  compression  ratio 

for  a  Chevy  5.74  liter  with  matched  power  output 
(early  intake  valve  closing) 
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Figure  14 
P-V  diagrams  of  Base  Cycle  and 
Modified  Cycle  having  matched 
power  output  (kW)  and  early 
intake  valve  closing,  2000  RPM 
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6.0  DISCUSSION  OF  IMPLEMENTATION  OF  A 
MORE-COMPLETE-EXPANSION  CYCLE 

The  reason  for  this  investigation  was  to  determine 
the  effects  of  the  more-complete-expansion  cycle  and  the 
feasibility  of  utilizing  such  a  cycle  in  the  application 
of  irrigation  water  pumping.   The  results  from  the 
present  study  appear  promising.   It  must  be  noted  that 
an  even  greater  advantage  can  be  achieved  in  an  engine 
where  the  displacement  volume  is  much  greater  than  the 
portion  used  for  output.   In  this  case  the  more-complete- 
expansion  cycle  will  not  only  use  the  entire  displacement 
more  efficiently,  but  can  also  cut  fuel  consumption  in 
half  while  reducing  power  output  by  75%.   The  same  is 
true  for  the  engine  used  in  this  study  if  it  were 
operated  at  half-throttle  or  below  rather  than  at  three- 
forths  throttle. 

Perhaps  the  question  remains  as  to  why  the  Chevy 
5.74  liter  engine  was  chosen  for  this  research.   As 
mentioned  before,  this  is  a  popular  engine  with  many 
modified  parts  available.   If  experimental  tests  are 
conducted  at  Kansas  State  University  in  the  future, 
this  engine  will  be  used  to  minimize  experimental 
costs.   Should  the  more-complete-expansion  cycle 
warrant  further  development  beyond  experimental 
testing  on  an  actual  Chevy  5.74  liter  engine,  a  more 
realistic  engine  can  be  purchased  and  modified. 
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6.1  Comparison  of  simulated  data  to  the  reference  case: 
Table  6  lists  the  percent  change  in  some  of  the 
important  performance  parameters  that  resulted  from 
the  engine  modification.   As  can  be  seen  in  the 
table,  all  of  the  modified  engines  exhibited 
approximately  the  same  power  output  as  the  reference 
engine.   The  modification  also  caused  a  13-20% 
increase  in  the  thermal  efficiency  which  was  expected 
due  to  the  increased  expansion/compression  ratio. 

The  results  shown  in  Table  6  are  necessary  to 
determine  the  effects  of  the  more-complete-expansion 
cycle.   Each  of  these  cycles  will  be  examined  in  the 
coming  sections  to  determine  the  cost  effectiveness 
of  implementing  the  modifications. 
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TABLE  6 

Comparison  of  the  modified  cases  to  the  reference 
case  for  determination  of  cost  effectiveness. 

Case  1:  matched  ign.  press.,  late  int.  closing 

Case  2:  matched  ign.  press.,  early  int.  closing 

Case  3:  matched  kW,  late  int.  closing 

Case  4:  matched  kW,  early  int.  closing 


Case  1 

Case  2 

Case  3 

Case  4 

% 

change 

in  kW: 

RPM: 

1600 

7.0 

9.4 

0.0 

0.0 

1800 

6.1 

8.3 

0.0 

0.0 

2000 

7.9 

7.0 

0.0 

0.0 

2200 

6.9 

6.1 

0.0 

0.0 

2400 

5.9 

4.6 

0.0 

0.0 

% 

change 

in  Eff: 

RPM: 

1600 

20.2 

13.  1 

17.5 

16.3 

1800 

19.4 

12.8 

17.1 

16.3 

2000 

19.4 

13.0 

17.0 

15.8 

2200 

18.5 

12.5 

15.7 

14.9 

2400 

18.2 

12.9 

15.7 

14.5 

% 

change 

in  FF  : 

RPM: 

1600 

-11.1 

-12.4 

-15.8 

-14.5 

1800 

-11.0 

-11.8 

-15.4 

-14.0 

2000 

-10.2 

-11.2 

-14.5 

-13.3 

2200 

-  9.7 

-10.8 

-14.  1 

-13.  1 

2400 

-10.2 

-10.4 

-13.6 

-12.6 
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6.2  Cost  analysis  of  engine  modification: 

An  important  feature  of  implementing  this  idea 
is  the  cost  of  the  modifications.   In  order  to 
incorporate  the  more-complete-expansion  cycle  on 
an  existing  irrigation  engine,  the  engine  must  be 
modified  in  the  following  ways.   The  compression 
ratio  needs  to  be  increased  to  the  appropriate 
value  and  the  cam  must  be  replaced  by  a  new  one 
ground  to  the  correct  specifications. 

The  compression  ratio  is  changed  by  planing 
the  cylinder  heads  down  to  reduce  the  volume  at  top 
dead  center.   A  price  estimate  of  $26  per  head  was 
received  from  Billy  Graham  Performance  Machine, 
Manhattan,  Kansas.   If  the  compression  ratio 
required  cannot  be  attained  by  planing  the 
cylinder  heads,  a  set  of  high  compression  pistons 
can  be  purchased  at  a  cost  of  $400. 

It  would  be  necessary  for  the  new  cam  to  be 
ground  by  a  professional  tool  and  part  shop.   A  few 
companies  were  contacted  to  determine  the  price  of 
the  cam.  (See  Appendix  II)   Each  of  the  shops 
estimated  that  a  new  master  would  be  needed  to  mill 
the  correct  profiles.   The  cost  of  this  procedure 
usually  runs  $340-350  per  profile  (i.e.  intake 
valving  and  exhaust  valving)  and  is  only  incurred 
once  if  several  cams  are  ordered.   If  the  company 
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happened  to  have  a  satisfactory  master  in  stock 
this  cost  would  be  eliminated.   The  other  costs 
include  the  cast  iron  core  at  $40-50  and  the 
actual  grinding  at  $250. 

Another  cost  created  by  the  modification  is 
the  replacement  of  the  gaskets.   The  gaskets  needed 
for  the  Chevy  5.74  liter  engine  are  priced  at  about 
$50.   Therefore,  the  complete  cost  of  the 
modification  with  the  two  required  masters  and  high 
compression  pistons  would  be  approximately  $1500. 
If  the  tool  shop  had  the  masters  in  stock,  the  cost 
would  be  reduced  to  $800.   Since  mechanic  rates 
vary,  a  generous  labor  cost  of  $800  will  be  assumed 
Therefore,  the  overall  cost  of  the  modification 
would  be  between  approximately  $1600-2300. 

Table  7  shows  the  amount  of  money  saved  per 
year  in  fuel  for  each  of  the  modified  cases.   An 
irrigation  engine  usually  operates  an  estimated 
1000-2000  hours  per  year.   If  the  engine  is  run  for 
2000  hours  per  year,  the  amount  of  money  saved 
would  be  significant  and  could  pay  for  the 
modification  in  less  than  a  year.   Whichever 
modified  engine  is  used  the  payback  period  would  be 
less  than  two  years. 

On  the  basis  of  this  analysis,  the  use  of  the 


TABLE  7 

Fuel  cost  saved  per  year  by  implementing  the 
more  complete  expansion  cycle.   Fuel  costs 
based  on  $0,097  per  Kg  of  methane  gas. 


RPM:  1600 
Hr  /yr 


operation 
1000 
1500 
2000 


Case  1 


854 
1281 
1709 


Case  2 


$  984 
1476 
1968 


Case  3 


$1220 
1831 
2440 


Case  4 


$1116 
1673 
2232 


RPM:  1800 
Hr/yr 


operation 
1000 
1500 
2000 


$  960 
1438 
1918 


$1043 
1565 
2086 


$1342 
2013 
2685 


$1220 
1831 
2440 


RPM:  2000 
Hr/yr 


operation 
1000 
1500 
2000 


$  994 
1491 
1987 


$1064 
1596 
2128 


$1412 
2118 
2825 


$1290 
1935 
2580 


RPM:  2200 
Hr/yr 


operation 
1000 
1500 
2000 


$1046 
1569 
2092 


$1100 
1650 
2200 


$1516 
2275 
3033 


$1412 
2119 
2825 


RPM:  2400 
Hr/yr 


operation 
1000 
1500 
2000 


$1202 
1804 
2404 


$1250 
1875 
2500 


$1603 
2404 
3206 


$1482 
2225 
2965 
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more-complete-expansion  cycle  is  warranted.   It  is 
also  recommended  that  this  idea  be  pursued  through 
the  experimental  testing  stage  to  determine  the 
actual  costs  and  advantages.   However,  different 
engines  must  be  examined  individually  by  a  computer 
simulation  or  by  experimentation  to  correctly 
incorporate  a  particular  engine's  features. 
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8.0  APPENDIX  I 

Engine  characteristic  parameters  for 
the  reference  Chevy  5.74  liter  engine. 


Fuel  Methane 

Bore ,  cm  10.16 

Stroke ,  cm  8.84 

Connecting  Link,  cm  14.50 

Compression  Ratio  8.50 

RPM  1600-2500 

Intake  Man.  Pressure,  kPa  76.00 

Intake  Man.  Temperature,  K  300.00 

Exhaust  Man.  Pressure,  kPa  101.40 

Lift,  cm  1.14 

Valve  Thickness,  cm  0.64 

Valve  Seat  Angle,  deg  45.00 
Valve  Diameter  (small) 

Intake,  cm  4.45 

Exhaust ,  cm  3.81 
Valve  Timing,  deg 

Intake  :  Open  (BTDC)  0.00 

Close  (ABDC)  76.00 

Exhaust:  Open  (BBDC)  77.00 

Close  (ATDC)  0.00 
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APPENDIX  II 

List  of  cam  grinding  companies  used  in  cost 
estimates  for  modifications. 


Ed  Iskenderian ,  Racing  Cams 
16020-T.S.  Broadway 
Gardena,  California 
(213)  770-0930 


Andrew  Products 
Chicago,  Illinois   6061 
(312)  992-4014 


E0NIC  Inc. 
466  East  Hollywood 
Detroit,  Michigan  48203 
(313)  893-8100 
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Appendix  III 

List  of  Variables  used  in  Computer 
Model  and  Complete  Program  Listing 
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LIST  OF  VARIABLES 


Units 


AA 

AEO 

AIO 

AREA 

AH(5) 

AL(5) 

B 

BB 

BH(5) 

BL(5) 

CD 

CM 

CP 

CPM 

CR 

CPR 

CPX 

CRIT 

CSUB 

CSUP 

CVR 

CUT 

CI 

CH(5) 

CL(5) 

DN 

DP 

DPH 

DPH1 

DPH2 

DT 

DV 

DW 

DWI 

DZ 

DTMAX 

DTMIN 

DH(5) 

DL(5) 

E 

EFF 

EK 

EO 

EVLD 

EVSD 

F 

FF 

FMEP 


fuel  equation  constant 
max.  flow  area  of  exhaust  valve 
max.  flow  area  of  intake  valve 
inst.  flow  area  for  valving 
gas  equation  constants 
gas  equation  constants 
bore 

fuel  equation  constant 
gas  equation  constants 
gas  equation  constants 
inst.  flow  coefficient  (drag) 
scaling  factor  for  chemistry  section 
inst.  constant  pressure  specific  heat 
constant  pressure  specific  heat  for  mixture 
compression  ratio 
reactant  gas  specific  heat 
residual  gas  specific  heat 
critical  pressure  ratio 
subsonic  flow  equation 
supersonic  flow  equation 
reactant  gas  specific  heat 
divisor  of  increment  equation 
constant  (2*L/s) 
gas  equation  constants 
gas  equation  constants 
fraction  of  burned  gas 
change  in  inst.  pressure 
change  in  pressure  due  to  heat  transfer 
pressure  change  of  hot  gases 
pressure  change  of  cold  gases 
change  in  inst.  time  increment 
change  in  inst.  volume  increment 
inst.  wotk  increment 
ideal  inst.  work  increment 
inst.  crank  angle  increment 
max.  time  allowed  for  integration  step 
min.  time  allowed  for  integration  step 
gas  equation  constants 
gas  equation  constants 
constant  (k-l/k) 
brake  thermal  efficiency 
constant  (1/3) 
energy  available  for  work 
exhaust  valve  large  diameter 
exhaust  valve  small  diameter 
fraction  of  do  loop  duration 
fuel  flow 

frictional  mean  effective  pressure 
92 


m 

m: 

m: 


cm 


kJ/kmol-K 
kJ/kmol-K 

kJ/kmol-K 
kJ/kmol-K 


kJ/kmol-K 


kPa 
kPa 
kPa 
kPa 
s 

3 

m 

kJ 

kJ 

rad 

s 

s 


% 

kJ 
cm 
cm 

kg/hr 
kPa 


LIST  OF  VARIABLES 


Units 


HC 

HC1 

HC2 

HRP 

I 

IVLD 

IVSD 

K 

KB 

KP 

KK 

KR 

L 

LFT 

LIFT 

M 

MDOT 

MEP 

MF 

MC 

MH 

MO 

MPS 

NC 

NCYL 

ND 

NM 

NMO 

NP 

NPO 

NW 

NX 

NXO 

N(l) 

M(2) 

N(3) 

N(4) 

N(5) 

P 

PC 

PE 

PH 

PI 

PM 

PNEW 

PO 

POLD 

POWER 

PR 

PX 

P1NEW 

PI 

P2 

P3 

P4 


inst.  heat  transfer  coefficient  kW/m'-K 

inst.  heat  transfer  coefficient  for  hot  gases  kW/m'-K 

inst.  heat  transfer  coefficient  for  cold  gases  kW/m'-K 

constant  pressure  heat  of  reaction  for  fuel  kj/kmol 

do  loop  variable  cm 

intake  valve  large  diameter  cm 
intake  valve  small  diameter 
specific  heats  ratio 

inst.  specific  heat  ratio  used  in  progressive  combustion 
specific  heat  ratio  for  product  gases 
do  loop  variable 

specific  heat  ratio  for  reactant  gases 

connecting  length  cm 

inst.  valve  lift  m 

max.  valve  lift  cm 

inst.  mass  in  cylinder  kmol 

change  of  mass  flow  w.r.t.  time  kmol/s 

mean  effective  pressure  kPa 

mass  of  fuel  kg 
number  of  carbon  atoms  in  fuel 
number  of  hydrogen  atoms  in  fuel 
number  of  oxygen  atoms  in  fuel 

mean  piston  speed  m/s 
do  loop  variable 
number  of  cylinders  in  engine 
do  loop  variable 

number  of  kmol  of  fuel  vapor  kmol 

number  of  kmol  of  fuel  at  stoich.  conditions  kmol 

number  of  kmol  of  product  gases  formed  by  comb.  kmol 
number  of  kmol  of  product  gases  at  stoich.  cond.  kmol 
counter  constant 

number  of  kmol  of  residual  gases  kmol 

old  number  of  kmol  of  residual  gases  kmol 

kmol  of  carbon  monoxide  kmol 

kmol  of  carbon  dioxide  kmol 

kmol  of  water  vapor  kmol 

kmol  of  oxygen  kmol 

kmol  of  nitrogen  kmol 

inst.  cylinder  pressure  kPa 

crank  case  pressure  kPa 

exhaust  manifold  pressure  kPa 

peak  cylinder  pressure  mPa 
constant  (3.14159) 

intake  manifold  pressure  kPa 

new  pressure  calculated  for  a  process  kPa 

atmospheric  pressure  kPa 

old  pressure  following  a  process  kPa 

engine  power  produced  from  cycle  kW 
pressure  ratio 

pressure  of  residual  gases  kPa 

new  pressure  before  compression  stroke  kPa 

pressure  before  compression  stroke  kPa 

pressure  after  compression  stroke  kPa 

pressure  following  combustion  kPa 

pressure  following  expansion  stroke  kPa 
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LIST  OF  VARIABLES 


Q  :  mass  flow  rate  through  valve 

QV  :  total  amount  of  energy  available  for  work 

R  :  universal  gas  constant  (8.314) 

RPM  :  engine  speed 

S  :  stroke  length 

SA  :  spark  advance 

SETANG  :  seat  angle  of  valve 

SFC  :  brake  specific  fuel  consumption 

SN  :  total  fraction  of  gases  burned  during  combustion 

T  :  inst.  cylinder  temperature 

TCOLD  :  inst.  temperature  of  cold  gases 

THNESS  :  valve  thickness 

THOT  :  inst.  temperature  of  hot  gases 

TM  :  intake  manifold  temperature 

TNEW  :  new  temperature  calculated  from  a  process 

TOLD  :  old  temperature  following  a  process 

TPOWER  :  total  power  output  for  all  cylinders 

TW  :  temperature  of  cylinder  walls 

TX  :  temperature  of  residual  gases 

T1NEW  :  new  temperature  before  compression  stroke 

Tl  :  temperature  before  compression  stroke 

T2  :  temperature  after  compression  stroke 

T3  :  temperature  after  combustion 

T4  :  temperature  after  expansion  stroke 

URP  :  constant  volume  heat  of  reaction  for  fuel 

V  :  inst.  cylinder  volume 

VBDC  :  cylinder  volume  at  bottom  dead  center 

VDISP  :  cylinder  displacement  volume 

VDOT  :  change  in  volume  w.r.t.  time 

VE  :  volumetric  efficiency 

VEO  :  volume  of  cylinder  when  exhaust  valve  opens 

VIC  :  volume  of  cylinder  when  intake  valve  closes 

VOL  :  old  volume  following  a  process 

WC  :  compression  work 

WCOMP  :  compression  work  done  during  combustion 

WC0MP1  :  compression  work  done  during  compression 

WEXP  :  expansion  work 

WL  :  work  done  during  valving 

WLOOP  :  work  done  during  valving 

WNET  :  net  work 

Y  :  air/fuel  ratio 

YCC  :  stoich.  air/fuel  ratio 

YMIN  :  min.  air/fuel  ratio  for  combustion 

Z  :  inst.  crank  angle 

ZC  :  combustion  duration 

ZCR  :  combustion  duration  in  radians 

ZDOT  :  change  in  crank  angle  w.r.t.  time 

ZE  :  exhaust  valve  opening  angle 

ZI  :  intake  valve  closing  angle 

ZIG  :  crank  angle  at  spark  ignition 

ZIO  :  intake  valve  opening  angle 


Units 

kPa/s 

KJ 

kJ/kmol-K 

rev-min 

cm 

deg 

deg 

g/kW-Hr 

K 

K 

cm 

K 

K 

K 

K 

lcW 

K 

K 

K 

K 

K 

K 

K 

kj/kmol 

m3 

m' 

m! 

m'/s 

% 

3 

m 

m5 
m1 

kJ 
kJ 
kJ 
kJ 
kJ 
kJ 
kJ 


rad 

deg 

rad 

rad/s 

deg 

deg 

deg 

deg 
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file  :  CrCLE 


uATFI  V 


KANSAS  3  T  -  T  c  UNIVERSITY  VM/i»  CMS 


JOB 


P=20fT=C2f)tNOLIST 


OTTO  CYCLE  (METHANE  =  U  E  L  ) 
AND  .REALISTIC  HALVING  FJR 


WITH  PC  J  ",  k  E  5  S  I  VE  COMBUSTION 
Thc  AN  IRsIoATI'JN  EUGINE 


SEAL  IV5 

REAL  N » N 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 


WCZ)=V3J 
G(Z)=2.« 

c  ( I )  =  V  D I 

acz;=pi* 


0,1 
M  ,  N 
,*, 
A,  A 
1  ,P 

eg, 

VLO 

TMA 

(5) 

'A 
G 

F 

A 
C-V 

SQR 

5P:c 


V  L  U  ,  K.  , 
P  ,  N  X  ,  N 
LtFiii 
0,c3,C 
2  >P*i  5 
A 1 0  $  C  U 
,  IVLUj 
x  ,CTMI 
,  AHC5) 
C  ZJ  =  V 
(Z)  =  A 

cz>  =  D 

CZ)  =  H 
OIiP*( 
TCC1** 

CSINC2 
2/<».+4 


i  <P 
.•IP 
S,T 
DD, 
T  M  , 
C  Vk 
FT, 
PIN 
LC5 
J  M  £ 
ITI 
I  VI 
T  T 
-Cu 
(SI 
2)/ 
WCZ 


,KP  , 
0,NX 
,  V,  Y 
DT,E 
Tw  ,  T 
,  N  M  0 
V3DC 
~A  ,  T 
}(BH 

AT 
0  N  A  L 
TlVc 
RANS 
SCZ) 
N(Z) 
GCZ) 
)/2 


LIFT,-,C,MEO 


, MH , MO , MPS 


<P, 
Tl, 
PM, 
OCT 
EW, 
)i3 
ANK 
ERM 
F  V 
R  A 
QRT 
*2) 
INC 


KR  ,  NM,h°  ,  N*',  NX  ,  PC  ,  ?  E  ,  ?I  ,  PM,  PO,PX 

T2  ,Th  ,  ZC, ZE , ZI 

YCC  ,  ZIO 

,ZJOT,WEXP 

VOISPtWLJOP.WCOMPl 

LC5),CHC5),CLC5),0HC5),DLC5) 

AW.GLE  Z 

IN  PCZ5 
0  L  U  K c  A c  T .  C  £  A  N  K  ANGLE 
k E  a  AT  C R A n <  ANGLE  Z 
CCl**2-CSINCZ))**2)-Cl)/2. 

2))/2. 


READ  IN  COEFFICIENTS  c  C  R  C-  A  S  PROPERTIES 

00  900  1=1,5 
*  0  0  R  c  A  D  (  5  ,  3  0  1  )  GmS,ALCI)i3LCI),CLCI),CLCI),A.iCI),:>HCI)»CHCI),DHCI) 
901  F0RMATCA4,=3.0,F6.2,=J.l,F3.2,c8.0,Fe.2,F9.1,F8.2) 

INPUT  OF  COEFFICIENTS  COMPLETE 

READING  J  -  F  u  =  L  UNO.  H  £  A  T I  \'  G  tf^U:  US 


IN  THc  CYCL: 


R  E  A  J  C  5  ,  1  7  5  )  G  A  6  ,  A  A  ,  ?  a  ,  M  C  ,  M  H  ,  M  0  ,  i  •<  ° 
175  FORMATCA6,F5.1,F7.*,3r*.0,=*.cn 


INPUT  0  i 


;UCL  DATA  COMPLETE 


FORMAT  jF  HEADING  fj^,  JUT  =  UT  OAT  A 

WRITE(5,202) 
Z0  2     FQR  m  ATC1H1  ,  2<  ,  1  rt  Y  ,  ix  ,  2HCR,  i\,  2MJA,  4X  ,  J^RP  n  ,  J  X  ,  ii  '*.  ;ii,  ,5X, 

C    12HVALVE  TI*!ilJG,<*X,3HMEP,5X,2'-,*W,7X,-.H3SFC,3X,7Hv0L  EFF, 

C     3X,2Ht=E,-DX,2H?;':,oX,2HPl,5X,2ki?rt,.*X,2rt.JL,*X,?nEF  =  / 

C    26X,6u'PISTjN,3X,3-'EXH,aX,JrtINT,5x,iJr<,PM,llX,9-"(;:-/Kv',.  -HR), 

C    HX,3Hi<>PA,5X,3H<PA,5X,5H^J4,^X,3nMCA,3X,2HKk,Fx,l'j;;/ 

C  26Xt5HS°E:0i*X,4ri0PEN»*X,5flCL0iE/ 

C  l2Xt7H0EGR5ES.7X,5HCM/5)»7X,7H0EGR&ES/ 

C  13X,-»H3T0C,l8X,iH3aDC,-X,^riAdDC) 

HEADING    COMB LETS 


PI  =  3  .141 5  9 
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file:  cycl 


w  a  t  f  :  v 


KJNSAi  iTaT:  UNIVERSITY  Vtf/SP  CMS 


R=e  .  314 


INITIALIZATUN    OF     i  N  G I N ' 


'ECIFICftTiaWS 


5=10.16/100. 

S  =  8. dn/1 0  0  . 

L  =  l4  .5/100. 

CR=3.5 

N  C  Y  L  =  3  . 

RPM=1S00. 

UPT  =  1  .14/100. 

TpiNc  SS  =  .  Sh/100  . 

SETANG=45.*PI/ldC. 

IVSD=4. ^5/100. 

EVSD=3. 61/100. 


ENGINE  SPECS.  CCMPLET 


C 1  =  2  .  * 

V  'J  1 3  P  = 
V50C=V 

V  T  0  C  =  v 
YCC=MC 

YMIN=Y 
URP=-H 
IVLD=I 

EVLD=E 

aiQsPi 

AEO=PI 


L/S 
PI* 

DIS 

BDC 
+  MH 
CC- 

P.P- 

vso 
vso 
/I. 
/I. 


B**2*S/4. 

P/C1.-1./CR) 

-V0I3P 

/  4  .  -  M  0  /  2  . 

MC/2. 

2430.*C«1C  +  Krl/2.-YCC-l  •  ) 

+2.*THNESS/T4NCSETANG) 

+2.*THNESS/TANtSET4NG) 

*(I  VLO+ 1 VSO  )*S'«'R  T  C  C  I VLD-1  VSO)** 2/4 


+  (LI-T-Th-j 


S)*«2) 


*  (  E  V  L  J  +  E  V  3  0  )  *  S  u  *  T  (  C  E  V  L  L  -  E  V  3  3  )  *  *  2  /  -  .  +  C  L  I  •=  T  -  T  H  N  E  3  S  )  *  *  2  ) 


iEGINMNG     U  =     CYCLE     CdLCULflTIGNS 


READ(5,400)     ZcD.ZiJ.ZlD 
^00    FORMATC3F10.0) 
DG     70  0     K  K  =  1 ,  1 0 
P  M  =  1  0 1 .  «♦ 
PC=250. 
?E=101.h 
PU=101.4 
T  M  =  3  0  0  . 
Ttf=423. 
RPM=RPM+100. 
ZDOT=2.*PI*R?M/60. 
Y=YCC+.2 
ZE=IEQ*PI/1B0. 
ZI=ZID*PI/180.+2.*PI 
VE0=W(-ZE) 
V  I  C  =  ,K  Z  I) 

P  1  =  P  M 

T1  =  TM 

NMQ=1.+4.76*YCC 

NP0=«C+MH/2.+3.76*YCC 

CPS  =  C4.76*YCC*2'i.3*AA  +  B3*Tl)/NMa 

CVR=CPR-R 
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CIL! 


CYCLE 


W  i  T  F  1  V 


«,  a  N  S  a  j  STAT;  UNIVERSITY  VVjo  CMS 


K,k  =  CPR/CVR 
14  N  M  =  P  1  •"  V  I  C  /  R  /  T  1 
VE  =  1  . 
N  X = 0  . 

:j  x  o  =  o . 

N  W  =  0  . 


20  N  W  =  N  W  + 1 

RN=NM/CNM+NX) 
CM=NM/NMO+NX/NPU 
IFCY-YCC)  3,2,2 

2  N(1)=0. 
NC2)=MC*CM 
N(3)=MH/2.*CM 
N(^)  =  (Y-YCC)=;=CM, 
MC5)=3.76#Y*CM 
GJ  TO  4 

3  N(1)=(2.*(YCC-Y))*CM 
NC2)=2.*CrCC-YMIN)*CM 
N(3)=MH/2.*CM 
N(4)=Q. 
N(5)=3.76#CM*Y 

4  QV=-CURP+(N(1)/CM)*281400.)*NM/N  1 C 
NP=NC1)+NC2)+UC3)+NC4)+MC5) 
ZC=40.+IS6.*(YCC/Y-1.1)**2+5.*CRPM/600.-1 

L0.'1?R6SSI  0  A    STROKE  3 r 3 INS 

Z=ZIQ*P 1/180. 
MCDMP1=0. 

CALL  COMP 

CjMPRESSICN  complete 


T=2000. 
6    TNEW=T-CUCT)-UCT2*CV/VT0Cja*CKR-l.))-^V)/CV(T) 
IrCAB3CT-TNEW)-5.)     *  »  *  ,  8 

8  T  =  T  'J  E  vi 
GO     TO     6 

9  T  3  =  T  N  E  U 

Pi=P2*(  tf/VTDC)**KR*NP*T3/(T2*I  V/VTOC  )**(KR-1  .  )  VCNM  +  NX) 

DELP  =  P3-P2*CV/VTQC)**K.R 

KP=(CVCT3)+R*NP)/CV(T3) 

PROGRESSIVE    C OM ? U S T I J N    B E  3  I  f>!  S    h E  k  E 

ZIG  =  Z 

ZCR=ZC*P 1/13  0. 

V=w( Z) 
P  =  ?2 

TC0L0=T2 

PQLD=P 

TOLO=TCOLO 
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ILc  :  CYCLE 


v,  A  T  F  I  V 


KANSAS  ST^T;  UNIVERSITY  V  M  /  S  ?  CMS 


>J  C  0  M  P  =  0  . 

W  £  X  P  =  0  . 

DPH1=0. 

0?H2  =  0  . 

DWI=0. 

N  C  =  <*  0 

?H  =  ?2 

DZ=ZC*PI/NC/130. 

0T=DZ/Zu3T 

S  N  =  0  . 

Du  300  1  =  1  ,NC 

Z=Z+DZ 

DV=W(Z)-V 

DN  =  PI*SIN((Z-ZIG)*PI/ZCRVNC/2. 

SN=SN+DN 

Ko=KR-CKR-KP)*SN 

D  P  =  (  -  K  S  *  P  -iO  V  +  0  E  L  P  -  V  T  D  C  #0  N  )  /  V 

DWI=DWI+CCP-PC)+DP/2.)*0V 

CALL  TEMPCQV,TC0L3,THQT,I,NC,DWI) 


HC1  = 
C  CW 
3PH1 
C?  =  C 
HC2  = 
DPH2 
DP  =  D 
0U  =  C 
IPC  J 
IFCQ 
P=P  + 
IFCP 
V  =  V  + 
TCCL 
PQLD 
TOLL) 
300  CONT 


303 


.0006 
ill)/ 

=  (P  +  3 
4.76=:= 
.0006 
=  (P  +  0 

P+OPH 

CP-PC 
w.LT. 

w.ST. 

OP 

.  G  T  .  P 

DV 

D  =  TOL 

=  P 

=  7C0L 

IIJUc 


*S*R?"I*(C  V(  THOT  )/NP+R  )v(P  +  jP)/THJT/^; 

V)v*C 1./2. ) 

P)=:«HClw(l.*I/NC)vACZ)*CT-A-THOT)*OT/CCV(THOT)«THOT) 

YCC*2  9.3+AA+cE*TCQL0)/Ntfa 

*  S*RPM*C?*(P+DP)/TCCLD/Rs:(W(ZI)/V)**(l./3.) 

P  )  *  H  C  Z  *  (  1  -  I  *  1  .  /  N  C  )  #  A  (  Z  )  *  C  T  vi  -  T  C  J  L  J  )  *  T  T  /  (  '  J  M  *  (  C  P  -  R  )  S  T  C  U  L  0  ) 

1  *  0  P  H  2 

)+D?/2. )*0V 

0)     WCOMP  =  WCOMP-0.i 

0  )    w  E  X  P  =  w  E  x  P  +  [■  w 

H)     ?h=P 

U*(P/PJLD)**C(KR-1.)/KR) 

D 


PROGRESSIVE  C0M3USTIJN  ClHPLcTc 

EXPANSION  STROKE  5 c  G 1 N $ 
CALL  EXP 

expansion  complete 
k  =  (cv(Th)  +  R:;=np)/Cv(T4; 

EXHAUST  ANU  INTAKE  3cGINS  -i  E  *  E 
CALL  EXHINT 

EXnAUST  AND  INTAKE  CCHPLfcTE 
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file:  cycle 


W  A  T  r  i  V 


KANSAS  STATE  uNIV-RSITY  V  M  /  3  P  CMS 


IFCA 

nxg  = 
pi  =  p 
r  1  =  r 

CPM  = 
CPX  = 
CPX  = 
CPR  = 
CVR  = 
KR  =  C 
GO  T 


SS(( 

NX 

1  N  E  W 

lNEw 

Ct.7 

NCI) 

CPX* 

C  C  P  M 

CPR- 

PR/C 

0  20 


NX0-NX)/NX).LT..01)  SO  TO  39 


6=PfCCM2  * .  i  +  4  A+5d*T  1  i*NM/NMQ 
*29.3+NC2)*37.2+N(3)v33.5+NC4)#29.3+N(5)*29.7 

NX/NP 

+  CPX  )/(NM  +  MX ) 

R 
Vk 


3  9  MF  =  CNM/NMO#(  MC*12.+MH*1  .+M0vl6.))*lU00. 
/c»NM#R#TM/PM/VOIS? 
SA=SA*130./PI 

FF  =  &0./2.*VIF*RPM 

MPS=2.*S*RPM/60. 

FrtEP  =  CC  .076*MPS#*2+.92*MPS+7. 235/14.  7  t  .  5  ),=  l  0  1 .  4 

i<NET  =  WEXP-WCOiMP  +  WLUGP  +  WCOMPl-CPHEPvVDISP) 

MEP=WNET/VOISP 

E Fe  =  -1 0 0  .  vWN  E  T*NMCI/ NM/ UR  P 

POWcR=rfNET*RPM/120. 

WL=WLQ0P*NCYL*RPM/120. 

WC=WC0MP1»NCYL*RPM/120. 

3FC=FF/POWtk 

T?QWER=NCVL*PUWeP 

CYCLc  CALCULATION  COM  =  L=.TE 

?H=?h/1000. 

PRINT  JUT  CYCLc  CALCULATIONS 

W  R I T  E ( 6 , 2  U  1 )  rfCR,SA,RPM,MPi»ZEO,2IO»MSPfTPCW  =  R,SrC,V=, 
C   P=  |PM, PI ,PH,lJLi ;FF 

201  FORMAT(F6.2,-5.1,F-3.1,F7.0,F7.2,2r?.l,F3.j,F:i.2,f:10.2,F6.2, 
C   3F8.1 ,F7.3.F5.  1,  =  7.25 


PRINTOUT  CCMPLzTc 


7U0  CONTINUE 
300  STOP 
END 


r'JNCTIQN 
INTERNAL 


U  (  X  )  USED 

ENERGY  0  c 


T  'J  :■  =  T  E  R  '4  i  N  : 
Th=  PROOUcT 


T*i 
C-ASE 


FUNCTION  U  C  X  ) 
REAL  N 
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I LE  :  CYCLE 


W  A  T  F I  V 


KANSAS  sTAT;  oNIVcRSITY  VM/5?  CMS 


COM  M  0  N 
C  0  M  M  0  N 
COMMON 
C  0  M  M  O  N 
COMMON 
COMMON 
COMMON 
U  =  0. 
IF(X  .G 
DO  2  I 
U  =  U  +  N  C 
RETURN 
00  <*  I 
U=U+NC 
RETURN 
END 


B,K,L»P»(*»R«3»T,V,Y,2 

AA|A0td8»Cl»0P»CTt£<»KS'»KRtN-M,NP,NW,NX,PC»?E»PIiPMlP0,PX 

P1,P2,P4,SA,TM,TW,TX,T1,T2,T-,ZC,ZE,ZI 

AEU,AlL'»CUT,CVR,NMO,RPM,YCC»tIO 

EVLD, IVLU, LIFT, V3DC,VQUT,Z00T,W EXP 

DTMAX,DTMIN,PlNEW,TlNEw,VQISP,WLGOP,WCaMPl 

NC5)fAHC&),ALC5),3HC3),ELC3)»CH(5),CLC5),0HC5),DLC5) 

T. 1600. )  GO  TO  3 
=  1  ,5 

I  )  *  C  A  L  C  I  )  +  (  5  L  (  I  )  -  R  )  *  X  +  C  L  (  I  )  v  A  L  DG  (  X  )  ) 

=  1.5 
I)*CAHCI)+C3HCI5-R)#X*CH(I)*ALOGCX); 


FUNCTI 

REAL  N 

COMMON 

COMMON 

COMMON 

COMMON 

COM  M G  N 

COMMON 

C  0  M  M  0  N 

U  R  =  J  . 

CP=(1. 

C  F  =  C  F  / 

UR=(CP 

R  t  T  U  R  N 

END 


FUNCTION  UR(X)  USED  TO  DETERMINE  THE 
INTERNAL  ENERGY  OF  THE  REACTiNT  GASES 


OfJ  urcx; 
M  ,  N  M  0 

3,K,L,P,*,R,S,T,V,ffZ 

AA,A0,s3|CltDP,0T,c'K,KP,KR,NM,N?,NW,NX,PC,PE,PI,PM,D3,PX 

P1,P2,P^,5A,TM,Tw,TX,T1,T2,Th,ZC,Z:,ZI 

AEQ,Aia,CUT,CVK,NMG,RPM,YCC,.:IG 

sVLD|IVLD,LIFT,V3CCfi/0aTtZ00Tt«,EXP 

07 M  AX, DTM IN, PI N  E  V* , Tl N  E  W ,VDIS?,WLOOP»  M COMF1 

NC5)fAri(i)tALC5)iBnCS>iSLC3).CrtC3)»CLC3)iDHC3)fDLC5) 

*  C  3  3  .  4  +  .  4  2  3  *  X  ;  +  Y  C  C  *  (  2  7  +  .  0  G  7  9  *  X  )  +  3  .  7  a  *  Y  C  C  #  C  2  7  .  b  +  .  0  0  5  1  *  X  )  ) 

NMO 

-  R  )  «  N  M 


FUNCTION  Cv(X)  US  EC  TO  FINj  T  H  ;  CONSTANT 
VOLUME  SPECIFIC  HEjT  OF  THE  PRODUCT  GASEs 


FUNC 
REAL 
COMM 
COM  M 
COMM 
COMM 
COMM 
COMM 
COMM 
CV=0 
IFCX 
00  2 


TION    CV(X) 

N 

ON  S,K,L,P,a,R,S,T,V,Y,Z 

Of,  J  A  ,  AD,  6  B  ,  CI  ,  OP,  OT  ,  EK,  KP  ,  KR  ,  NM,NP,N*  ,NX,  PC,Pr  ,  PI  ,  FV  i?0  .  ?X 

ON  PI,  P2,  F4,  S  A,  TM,Tw,TX,Tl,T  2,  T-*,ZC,Zi,Z  I 

ON  AEO,Aia,CUT,CVR,NMQ,RPM,rCC,ZI0 

ON  EVLD,IVLD,LIFT,VbOC,VDOT,ZUOT,rtExP 

OfJ  D TMAX  ,  0  TM  IN  ,  P  IN c W  ,  T 1  N  E  W  ,  V 0 1 S  P  ,  WLUC?  ,  w C 0 MP  1 

ON  NC53,AHC5),ALC5),EHC5),2L(5),CriC3),CLC5),0riC5),0LC5) 


.GT  .  l&OO  .) 
1=1,5 


,0     Ti 
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=  lLfc  :  CYCLE 


W  A  T  F  I  V 


KANSAS  STATE  UNIVERSITY  V  M  /  S  P  CMS 


10 


cv  =  cv  +  nci)=::(Bl(I)-r  +  cl(I)/x) 

RETURN 

0  0  <♦  1  =  1,5 

CV=CV+NCI)*CEHCI)-R+CH(I)/XJ 

RETURN 
END 


SUBROUTINE  TEMP  USEO  TO  FIND  ThE  TEMPERATURE 
OF  THE  HOT  ftf-jC  ClLO  SASES  DURING  TriE 
PROGRESSIVE  COMBUSTION  PROCESS 


SUBR 
RcAL 
COMM 
CO, MM 
CO  MM 
COMM 
CUMM 
COMM 
CGMM 
F  =  l. 
TCJL 
501  = 
E0  =  ( 
TT  =  2 
ThEA 
IF(A 
TT  =  T 
GU  T 
TmCT 
RETU 
ENO 


OUTI 
K  ,  K 
ON  5 
ON  A 
ON 
GN 
ON 
ON 
ON 


91  /N 

D  =  TC 
UCTC 
EQ1  + 
000. 
=  TT- 
BSCT 
N£w 
0  6 
=  TNt 
RN 


NE  TEM=(Qv,TCOLl;,THl.iT,I,NC,OWI) 

P,KR,M,MQaT,NM0,NM,NP,NX 

,K»L,P,GtR»S,T,V, Y,I 

A,A!j»33,Cl,GP,0T»=<,KP,K.R,NM,NP,NW,NX,PC,FE,PI,PM,Pd,PX 

1,P2,F*,SA,TM,TW,TX,T1,T2,T-,ZC,ZE,ZI 

EC,AlU,CUT,CVR,NMO,KPM,YCC,Zia 

VLD,IVLD,LIFT,  VBPC , VDQT , ZUQT , WEXP 

T  M  A  X  ,  l  T  M  I  N  ,  P  1  N  E  W  ,  T  1  N  E  W  ,  v  CI  S  P  ,  W  L  J  0  P  ,  W  C  0  M  P  1 

C5),AHC5),ALC5),3HC3),3LC3),CHC5),CLC5),DH(b),0LC5) 

C 

0  L  C ■  -::  (  C  P  + : IP  )  /  - )  ■■■  *  C  C  <  K  - 1 .  )  /  n  R  ) 

0  L  D  )  +  J  V 

F*UR  (T2;::=T2+=->;v-F:::u(TCCLu)-F:,:.;v'-C1UI) 

(UCTT)-;D J/CVCTT) 
T-TNEw)-5.)  9,9,3 


A 


SUBROUTINE  CQ-MP  US  EC  TO  SIMULATE  THE 
COMPRESSION  PROCfcSS  WITH  HEAT  TRANSFj 


SUER 
.REAL 
COMM 
COMM 
COMM 
COMM 
CGMM 
COMM 
CGMM 
W(Z) 
A(Z) 
NC  =  5 
SA  =  C 

oz  =  c 

0T  =  0 


0  U  T I N 


ON 
ON 
ON 
ON 
ON 
ON 
ON 


b  , 
AA 
PI 
AE 
EV 
OT 
N( 


=  V  B  U  C 
=  P I  *  B 
0 
ZC/2. 

PL-l- 
Z/ZOO 


C  0  H  = 
R  ,  M , 
L,P, 
U,  LB 
2,F4 
Alu  , 
iIVL 
X  ,GT 
,AHC 
OISP 
2/4. 


M  LOT,  N  M  ,  N  M  0  ,  N  X 
w,R,S,T,tf,Y,Z 


,C1 


T , iK.KP  ,  K?,riM,f!F  ,  N  W  ,  N  X  ,  ?C  ,  P  E  ,  =>  1  ♦  P  M  ,  °  0  ,  P  X 


,SA,TM,TW,TX,Tl,T2,T«t,ZC,ZE.ZI 

CUT,CVR,NKO,'K?M,YCC,ZIO 

0, LIFT, ViDC,V OUT, ZDCT.W EXP 

MI N , P IN  =  w , T 1NE  M , V  0 1 S  P , W  LOTP , WCU  MP 1 

53,ALC3),BHCS),3LC5),CHC5),CLC5),0HC5),0LCS) 

*Cl.-CQSCZ)+SQRTCCl**2-CSINCZ);**2)-Cl)/2. 

+4.*W(Z)/B 


)*PI/1S0. 
SA)/NC 

T 
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file:  c  r  c  l  e 


WATFI  V 


KANSAS  5  T  A  T  d  'J  K  I  V  b  R  S  I  T  V  VM/SP  CMS 


P  N  E  W 

T  N  e  w 

VCL  = 
DPH  = 
00  1 
Z  =  Z  + 
V  =  W( 
P  =  PN 
T  =  P* 
CP  =  C 
HC  =  . 
0PH  = 
DV  =  V 
0P  =  C 

ow=c 

WCOM 

p  n  e  w 

VGL  = 
TNEW 
CONT 
P2  =  P 
T2=T 
RETU 
END 


=  P1 
=  7  1 
W(Z) 
0. 

1  =  1 
OZ 
Z) 

E  W  *  C 
V/NM 
<*.  7  b 
0006 

P:.=  HC 

-VOL 

P-PN 

(PNE 

P1=W 

=  PNE 

V 

=  PNE 

INUE 

NEW 

NEW 

RN 


» i  J  C 


VGL/V)**KR 

/R 

*  r  *  2  9  .  3  ♦  A  A  ♦  E  3 :;:  T  )  /  N  M  G 

*S*RPM*CP*P/T/R*CW(ZI)/V)**(l./3.) 

SA(Z)*(TW-T)*DT/(NM*CCP-R)*T) 

=  W  )  +  OPH 

W-PO  +  OP/2.  )*0V 

C  <J  M  P  1  +  0  w 

W  +  OP 

W*VOL/R/NM 


SUBROUTINE  EXP  U3E0  TO  SIMULATE  THE 
EXPANSION  PROCESS  wITH  w  =  a  y  T=>ANSF  =  R 


sues 

REAL 
COMM 

COMM 
CO  MM 
CGMM 
COMM 
COMM 
COMM 
W(Z) 
A(Z) 
N0  =  o 
OZ  =  C 
OT  =  D 
VCL  = 
OPH  = 

z  =  z  + 

v  =  wc 

P  N  E  W 
TNEW 
HC  =  . 
DPH  = 
DV  =  V 
DP  =  C 
0  w  =  C 


OUT 
Kt 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

=  V3 

=  PI 

0 

2.=:= 

Z/Z 

V 

0. 

OZ 

Z) 

=p* 

=  PN 
000 
PNE 
-VO 
PNE 
CP- 


INE 
KP. 
d,  K 
^, 
Pit 
AEO 
EVL 
DTM 
IK  5 
OC- 
*  3  * 


EXP 
K  k  »  M 
.L,P 

Ad,  b 
32,r 

» A  I  0 
0,  IV 
AX  ,  0 
),  Art 
VOIS 
*2/«* 


,  MOOT,  NM  ,NP  ,'U 

,Q»RiS»TfV,r,Z 

t  ,  C  1  ,  0  P  ,  0  T  ,  E  <  ,  K,  P  ,  i\  R  ,  N  ,-1 ,  N  p  ,  M  .V , ,  j  <  ,  o  c  ,  ?  E  ,  P  I  ,  P  M  ,  3  J  ,  P  X 

A,Sa,TM,T.v,TX,71,72,"!"H,ZC,Z-,ZI 

»  CUT  » C  V  * , NMQ , R PM , YCC  >  £  I C 

LDtLIrT,V6DCiV0UT,Z'JCT,WEXP 

TM I H  t  P  1 N  5  W  »  T 1 N c  W  »  V  0 1  S P  »  WLU 0 P  i  W C Q  w P 1 

CJ)fALC53f3rtC5)»3LCD)fCriC53,CLC53iOHC5).DL<5) 

P* C  1. -C GS  (  Z  )  +  S3RT(C  1**2- (S  I, \(Z  ))**?.) -Cl)/2. 

.+4.*W(Z)/B 


PI-Z- 

OOT 


Z  c  )  /  N  , 


CVOL/V)**CKP) 
EWvV/NP/R 

6  *  S  *  R  P  M  *  (  C  V  C  T  N  c  W  )  /  N  P  ♦  R  )  *  P  N  c  W  /  T  N  E  W  /  R  *  (  W  (  Z  I  )  /  V  )  **  ( 

W *H C*  A  C  Z  ) *  C  T  'W - T N E  W  )  ~.  D  T /  C  C  V  (  T M  E  *  )  *  T  N  5  A  ) 

L 

W-P >+DP H 

PO  +  CP/2.  )*OV 


1./3. ) 
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FILE  :     CYCLE 


ft  a  t  p  i  v 


KANSAj     STAT:     UNIVERSITY     V-1/SD     CIS 


WEXP=WEXP+Qfc 

P=P+OP 

VOL  =  V 

T  =  P;.:\/uL/K/NP 

IF((Z+DZ).LT.(2*PI-ZS))  G3  TO  1 

T<f=T 

P4  =  P 

RETURN 

EfJD 


SUSfiiJUTINS  EXHINT  USED  AS  THc  CONTROL 
STRUCTURE  FG«  THE  VALVING  ROUTINES 


SU5R 
CQMM 
COMM 
COMM 
COMM 
COMM 
COMM 
CCMM 
CUT  = 
OTMA 
OTMI 
WLOO 
EK  =  . 
l  =  -l 
CALL 
CALL 
CALL 
RETU 
END 


iUTI 


ON 
ON 

0  N 

ON 

CN 

ON 

ON 

10. 

X  =  3. 

H=.c 

P  =  0. 

333 


N  E  EXHINT 

,K tLtP»CiRiS»T,V,Y  iZ 

ft , AO , b3  ,  CI , D  P »  C T  »  =  < , KP  » K  R  »  MM , N P , N W  »  NX , PC , F  E , PI t  PM  f  P J  ,  P X 

l,P2»P«nSApTM,Tw,TXfTl»T2,T'*,IC,Z  =  ,ZI 

EOiAIO»CUT,CVRfNMO,RPM,YCC»ZID 

VLD , I V LD  f  LI *T , V3 DC  t  VDOT , Z  JOT , WE X  P 

TMAXfOTMIN,PlNEW,TlNEW,vOISP»WLUOPt^COMPl 

CS)tAHC5)fALC5)ibriC3)»flLC3).CriC5).CLC53tDriC5),0LC5) 

'•PI/130. /ZDUT 
5=;=  PI/180. /ZUOT 


EXH 
TDC 

sue 

RN 


SUaRUUTINE  EXr\    USErJ  T  rJ  SiM'JLATc  T*i 
EXHAUST  PROCESS  DURING  ThE  CYCLE 


SU3R0U 
REAL  K 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
CUMMON 

gcz;=2 

F  C  I )  -  V 

-KZ)=V 
c=CK-l 

CRIT=( 
CSUd=- 


Tint   t 

,LFT,L 
i,K,L 
AA,  AO 
PI  ,P2 
AEO  »A 
EVLO, 
U  T  M  A  X 
N  C  5  )  , 
.  *  S  Q  R  T 
DISP«C 

eoc-vd 
.  )/«, 

U  +  1.3 

5wRT(2 


Xh 

I  F  T  ,  M  , 

»  P  ♦  Q  t  R 

iEB.CI 

»P4,Sa 

IQ.CJT 

IVLC  ,L 

i  JTHIiJ 

AHCrO, 

CC1**2 

SINC2. 

ISPvCl 


MDOT 

i3  iT 

lOPt 
•  TM, 

,  CVR 
IFT, 

,pk; 

AL(5 
-(SI 
*Z)/ 

.-co 


,  N  M  ,  N 
»  V  ,  1 , 
D  T  ,  c  K 

T  W  ,  T  X 
» N  M  3  , 

V30C, 

e  w » r  i 

) ,  3  ri  ( 

NCZ)) 

GCZ)- 
SCZ)  + 


,<P,KR,NM,NP,sJ^,NX,PCi?F,PI,PM,PO,PX 

iTl»T2»T»iZCtZEf ZI 

R3M,tCC,ZIG 

V90T,ZC0T,W=XP 

NEW.VDISPfWLGOP.WCGMPl 

5)f3LC5),CH(5),CLC5),DM(D),OLC3) 

**2) 

sincz)  )/:. 

o 1 R  T  C  C 1 ** 2-C  5  I  IK  Z ) ) *w2 )-C 1 ) /2 . 


/Z.  }=s*C</CK-l. )) 

0O0./E/R/J3. 97) 
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file  :  CYCi_e 


W ATFI  v 


KANSAS  STAT:  UNIVERSITY  VM/SP  C  .-1  S 


C5UP 
CSU? 
?  =  P<+ 
T  =  T* 

v  =  «i  c 
M  =  P  $ 

AREA 
LFT  = 
C0  =  - 
VOOT 
IFCP 
DP  =  C 
IFCP 
PR  =  C 
MOOT 
GO  T 
IFCP 
PR  =  C 
MOOT 
GO  T 
MOOT 
OP  =  C 
0  =  K* 
GO  T 
MOOT 
Q  =  K* 
CALL 
P  =  P  + 
Z=Z  + 
V  =  WC 
T  =  T4 
M  =  P--;: 
IFCZ 
?X  =  P 
TX  =  T 
NX  =  M 
RETU 
ENO 


=  C  2  .  /  C  K  + 1 .  )  )  **(  C  K  + 1 .  )  /  C  K. - 1 .  )  ) 

=-5gRTCCSUP*<*lU00./2S.*7/R) 


Z) 

V/T/R 

=AE0*(A3SCSINCPI*(Z+Zc)/(PI+ZE)))) ** E  K 

L I  F  T  *  (  A  6  S  (  S I  M  (  P  I  *  (  Z  +  Z  E  )  /  C  P  I  -  Z  c  )  )  )  ;  *  *  E  \ 

20./3.*CLFT/£VLO-.3)**2  +  l./3.-CLPT/!EVLO-.3)  +  .7 

=FCZ)vZOOT 

•LT.Pc)  GO  TO  4 

P-PE5/CUT 

/PE.GT.CRIT)  GO  TO  2 

P/PE)**c 

=CD*AREA*CSU3*SjRTCPRvCPR-l.)/T)*PF 

0  3 

E/P. GT.CRIT)  Gj  TG  5 

PE/P)#«E 

=-CD*AR=4$CSU3*SQRTCPR*CPR-l.)/T)*PE 

0  6 

=-CD*AR5A*CSUP*P5/SuRT(T) 

PE-P)/CUT 

CT*MOCT*R-PwVOOT)/V 

0  7 

=  C  0  *  A  R  E  A  *  C  S  U  P  *  p  /  S  w  R  T  (  T  ) 

P*(M0OT/M-V3OT/V) 

STEP 
OP 

Z  C  0  T  *  D  T 
Z) 

*(P/P  ■+)**£ 
V/T/R 
.LT.PI)  GC  TO  1 


RN 


SUBROUTINE  T  J  C  USE  J  TO  D  E  C  k  E  u  S  z  THE 
CYLINOER  PRESSURE  3  E  L  0  W  THE  MANIFOLD 
PRESSURE  SO  THwT  ThE  I  N  TAKE  STROKE 
CAN  6  E  G I N 


SUBROUTINE  TOL 

REAL  IVLO,K,LrT,LIsT,M,MDOT,Mx 

COMMON  3,K,LtPiitRiS»T,V,Y,2 

COMMON  AA,A0,i5,Cl,0Df0T»EK,KPi!<R»NM,NP,NW,i>IX,PC»PE|PItPM,PJiPX 

COMMON  Pl,P2tP4,SA,TM,TW»TX,Ti,T2,TH,IC,Z=iZI 

COMMON  A  E  U  ,  A  I  G  »  C  U  T  ,  C  V  R  ,  N  M  0  ,  k  P  ;■',  ,  T  C  w  ,  Z  I  j 

COMMON     EVLD,IVLD,LIFT,V3  0CiV'JGT,ZOClTtWcXP 

COMMON  DTMAX  ,OTMIN,  PlNFJ  ,TlNE*J>  VCISP  ,  WLGOP,WCO""=>l 

COMMON  N(5),AH(5),AL(5),8H(b),ELCa),CM(5),CLC5),OHCi),OL(!5) 
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FILI 


CrCLE 


USTFIV 


K  a  N  5  A  S  S  T  A  T  •;  UNIVERSITY  V  ;1  /  S  P  C  -1  S 


GC 

FC 
rfC 

IF 
IF 
RE 
*  10 
NC 
V0 
OZ 
DJ 
Z  = 

V  = 
?N 
DP 
DV 
0*1 
WL 
VO 
P  = 

5  CD 
RE 

1  £  = 
CS 
T  = 
P  = 
M  = 

V  = 

2  DP 
PR 
AR 
LF 
CO 
MD 
VO 
Q  = 
CA 
P  = 
Z  = 

V  = 
T  = 
M  = 
IF 
Pa 

CS 

3  OP 
PR 
AR 
LF 
CO 
MO 
VO 
3  = 
Cm 
P  = 


=  2.-5 

aVOI  S 

=  V3DC 
1 0  .  N  E 

M  .  L  T  . 
RN 

I-t-ZIO 
0 

wcz) 

ZU-PI 

1*1. 

OZ 
Z) 

=  P  *  c  v 

NEW-P 

-VOL 

(P-PC 

P  =  WLO 

V 

DP 

INUE 

RN 

-1.)/ 

=  -SJR 


JRTCClv*2-CSINCZ))**2) 
P*(SiNC2.*Z)/G(Z)-SINCZ))/; 
- VD IS P*( 1 . -COSC  Z )* SSR T ( C 1*: 

.  0  .  )  ;  (j  T  D  4 

P)  Go  TO  1 

*  P I  / 1  5  0  . 


)  /N  C 

NC 


2-CSINCZ))**2)-Cl )/: 


OL/V  )«:.:< 


)+0P/2)*3V 
OP+OVn 


K 
TC2 


Q0./E/23. 97/S) 


P  M  )  /  C  U  T 

PM)**S 

1 0*C  AdS C  3  IN ( P I *C  Z- Z I ) /C  P I- Z I ) )  )  ) #*  EK 

FT#C  ABS C  S IN C P I *( Z-Z I )/ C  P I-Z 1 ) 3 ) )**EK 

./3.*CLFT/lVLU-.3)ft*2+l./3.*CL=T/IVL0-.3>- 

D*AREA*CSU9"S=jRTCPRvCP^-l.)/T)«? 

=F(Z)*ZDQT 

P*(MOUT/M-VOCT/V ) 
STEP 

DP 

ZDQTssDT 

Z) 

*< 

v/ 
-p 


P/PX 
T/R 

M.GT 
fc3 

csua 
-p)/ 

/?):.= 

I  0  *  c 
F  T  -,'•-  ( 
./3. 

0  -.=  A  K 

(  Z  )  y 
M*MD 
TEP 


.0.63)     G  G     TO 


CUT 


A  3  S  (  S  I N  C  P  I  *  C  Z-  Z  I )  /  C  P  I  -  Z  I  )  )  )  )  *  v.  E  < 

AdSCSINCPI*CZ-ZI )/CPI-ZI))))**ER 

$CLFT/IVLC-.3)**2+l./3.«(LFT/IVL0-.3)+.7 

E  A  *C  S  U  3  :.'•  5  <<  R  T  (  P  5  *  (  =  R  -  1 .  )  /  T  M  )  *  t>  -1 

Z3QT 

QT*R-P*VJUT)/V 
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ILc:  CYCLE 


w  a  t  i= :  v 


*  -  N  S  A  i  S  T  J  T  £  UNIVERSITY  VM/SP  CMS 


z  =  z  + 

V  =  W( 

N|  =  M  + 
IFCP 

IPC? 
RcTU 
DP  =  < 
Pft  =  C 
AREA 
LFT  = 
C0  =  - 
T  =  TX 
MOOT 
VOOT 
Q=K# 
CALL 
P  =  P  + 
Z  =  Z  + 

v=wc 

M  =  M  + 
IFCM 

RETU 
END 


2  D  0  T  * 
I) 

MOOT* 
.LT.P 
.bT.P 
RN 

P-PM) 
P/PM) 
=  AIQ* 
L  I  F  T  * 
20. /3 
*  C  P  /  ? 
=  -CD* 
=  FCZ) 
P*(MO 
STEP 
DP 

ZDOT* 
Z) 

MOOT* 
.LT.N 
RN 


DT 


DT 


M  .  A  M  D  .  M  .  L  T  .  N  X  )  G  J 
M.ANC.M.LT.NX)  G  'J 

/CUT 


TO 

TO 


(AaS(SIN(PI*(  Z-ZI  V(PI-ZI)  )  ))**£* 

C4  33C3INCPI*CZ-ZI)/CPI-Zi))))*v£K 

.*(LFT/IVLD-.3)**2+l./3.*CLFT/IVL0-.3)+.7 

X)**E 

AR E A*C SU3*S3RT( PR* ( PR -l.)/T)*o 

vZJOT 

3T/M-VD0T/V) 


UT 


DT 

X)  GQ  TO  6 


SUBROUTINE  SJC  USED  TO  SIMULATE  THE 
INTAKE  STROKE  DURING  THE  CYClE 


SUBROUTINE  SUC 

REAL  I  V  L  D  ,  K  ,  K  K  ,  L  F  T  ,  L  I  F  T  ,  M  ,  '-1  0  0  T  ,  N  '-1  ,  N  X 

COMMON  B,k.»L»P»2»R.S»T,V,Y,2 

COMMON  4A,AU»BefClfDP|CT»SKtKP»KR,NMiNPfNW»NX,PC»PE,PI,PMfPJ,PX 

COMMON  Pl,?Z|P4,SA,TM,Tw,TX,TifT2tT<*,ZC»Z5»ZI 

COMMON  AE0»AlOiCUT,CVR,NMG,R?M,YCC,Zia 

COM  M  ON  EVLD,lVL0fLlFT,Ve0C,V0aT,ZC0T»'W5xP 

COMMON  DTMAX , DTMIN  ,  P  1  N E W  ,  T 1 N E W  ,  V 0 i  S D  ,  -JLOOP  »WCUMF 1 

COMMON  NC5),AHC5),ALCp)»3H(5)f5LC5)»CHC5),CLC5),DH(5),DLC5) 

G  C  Z ) =2 . *S  QKT  C  C l**2-( S I N( Z ) )#*Z ) 

FCZ)=VDISP*(SIM(2.*Z)/GCZ)-5IN(Z))/2. 

W(Z)=VSUC-VDISP*C1.-C0S(Z)+S3RTCC1**2-CSI;KZ))**2)-C1)/I. 

M=NX 

K=Kk 

E=(K-1. )/K 

CRIT=C(K+l.)/2.)  ** C  K / ( K -  1 . ) ) 

C SUP  =  C  2  .  /  C  K  + 1 .  )  )S=?C  <  K  + 1 .  )/  CK- 1  .  )  ) 

CSUP=$SRT(CSUP*K*1000./23.97/R) 

CSU6  =  oCRT(.2000./E/2o.97/R) 

V=UCZ) 

L  =  l 

AREA  =  AIC*CA3SCSINCI*CZ-ZI)/CPI-ZI))))**EK 

LFT=LIFT*CA3SCSIMCPI*CZ-ZI)/CPI-ZI))))**=N 

CO=-2U./3.*CLFT/IVLO-.3)**2+l./3.*(LFT/IVLJ-.3)+.7 

VDOT=F(I)*ZOOT 

GO  TO  C  2  ,  3  )  t  L 
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pile:  cycle 


A  A  T  r  I  V 


TA' 


UNIVERSITY  V  M / 5  P  CMS 


2  IFCZ.LE.2.*PI)  oO  !u  J 
IFCP.LT.PM-.6ti)  JG  TO  4 
L  =  Z 

P  =  P  M  +  .  6  3 

3  IF(P.LT.PM)  G  Q  TO  fc 
OP  =  (  p-P,v|)/CUT 

I=<P/PM.GT.CRIT)  GO  TO  d 
PR=(P/PM)**E 

MDaT=-CU*MRE4s.5CSU8*SQRTCPR*CPR-l.J/T)*P 
GO  TO  9 
3  MDOTa-CDSAREASCSUP^P/SQRTCT) 

3  QsK*P#CMDQT/M-Vl)OT/V) 
GO  TO  7 

4  DP=(PM-P)/CUT 
IFCPM/P.GT.CRin  GO  TO  5 
PR=C PM/P)**E 

MOOT  =  CD*AREA3CSU3*S«,RT(PR*(PR-l.)/TM)*PM 
GO  TO  6 

5  MOOT  =  CD*AR£A*C5UP*!,M/S'3RT(TM) 

6  Q  =  K* C T M*MU DT*R-P* VOQ T ) / V 

7  CALL  STEP 
P  =  P  +  0  P 
Z=Z+ZUOT*UT 
V=W(Z) 
M=M+MDOT*DT 
IF(Z.LT.ZI)  GO  TC  1 
TlNcW=P*V/R/M 

P  1  N  E  w  =  ? 
N M  =  M -NX 
RETURN 
END 


SUBROUTINE  STEP  USED  T  Z    DETERMINE  TH: 
APPROPRIATE  STcP  SIZc  FOR  THE  INTEGRATION 
TECHNIQUE  IN  TrE  VALVING  SU*  -C  GU  T  I  Ns  S 


SU3R 
C  0  M  M 
COMM 

COMM 
COM  M 
COMM 

COMM 
C  0  M  A 
IFCQ 
QT=- 
IFCD 
DT  =  0 

DP=y 

GO  T 
IFCD 
OP  =  - 


OUTI 

ON    J 


ON  A 
ON  P 
ON  A 
ON  E 
ON  U 
CfJ    N 


N  E     S  T  =  P 

i  H  i  L  »  P  i  ,  ,  ?  ,  5  i  T  i  V  i  y .  I 

A,ADt  =  3,Cli0P,0TfEKi<PtKR,N*!,NP,Nw,.MX,PC»PE|PI.PM,P0»PX 

lfP2iP*iSAtTM,TW,TXfTl,T2,T<»»i:CtZE.II 

EOiAlOtCUT,CVR|NMOtRPM,YCC,ZI'J 

VL0,IVLJ,LIFT,V30C,V!:GT,Z0GT,^EXP 

T  M  A  X  ,  J  T  M  1  N  ,  P  1  N  E  w  ,  T  1  N  E  w  ,  v  0  I  S  s  .  w  L  C  0  F  ,  w  C  J  V  ?  1 

C33fAHC5)i4LC5)fSriC5>fSLC5)fChC5>»CLCS}|0riC5)i0LC5) 

5  ,  6 


)  2, 

CP/J 

T  .  L  T  .  0  T  M  A  X  )  GO  TO  ? 

TMAX 

*DT 

0  9 

T.LT.DTMIN)  G  J  TO  4 

DP 
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file:    cycle 


W  A  T  F 1 V 


jjj     S  T  x  T  t     UNlViRilTY     V  M  /  3  P     CMS 


SJ     TO     9 

DT=OTMIN 

D  P  =  J  *  D  T 

GO    TO     9 

0?  =  0. 

DT=DTMAX 

GO    TO     9 

OT=OP/S 

IFCDT.LT.OT.IAX)     GO     TO     c 

D  T  =  D  T  M  A  X 

DP=U*OT 

GO    TO    9 

DT'DTMIrt 

DP  =  il*DT 

GO    TO     9 

IFCDT.LT.t5TM  IN) 

WLOOP=WLOQP+CCP 

RETURN 

END 


GO     TO     7 
PO+DP/2.  )WDO 


JENTRY 
CO     299180 
C02     5683= 


H20 

02 

N2 

CH* 

77. 


d8  923. 

4338o. 

31317. 

20.1 


37.8? 
66.  1  7 
49.3*- 
42.27 
37. ->6 
0.0051 
0. 


-4571. 9 

-1163*. 0 

-7:>40.  3 

-6635. <+ 

-4559.3 

1.       <*. 

76. 


-31.10     3 
-200.00 
-117.00     1 
-i>5.l5    1 
-34.82 
0.  797 


T*OT 


0  9  0  7  0 

5-670 
2  7  0  10 

4  <♦  e  5  9 

5  70 


39.29 
ts».  58 
=  0.4j 
46.25 
3  3.32 


-6  20 
-1  =  97 
-1921 
-lo79 

-o75 


1  .9 

9.0 
2.0 
3.0 
*  .  4 


4  2.77 
20.^0 
0  4  . «-  0 
92.15 
50.24 
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ABSTRACT 


The  more-complete-expansion  cycle  can  be  applied 
to  a  stationary  irrigation  engine  commonly  operated  at 
constant  partial-load  burning  liquid  fuels.   The 
analysis  involved  the  development  of  a  digital 
computer  model  to  simulate  a  real  engine. 

The  Chevrolet  5.74  liter  (350  cu.  in.)  engine  was 
simulated  to  determine  the  performance  changes  and 
possible  advantages  of  applying  the  more-complete- 
expansion  cycle  features.   The  fuel  used  in  this  model 
was  methane  gas.   The  generated  performance  data  is 
presented  in  figures.   The  modified  engine  specifications 
are  results  of  these  figures  and  are  used  to  determine 
the  cost  effectiveness  of  the  concept. 

The  performance  results  in  conjunction  with  a 
cost  analysis  suggest   that  a  fuel  saving  greater  than 
11%  can  be  achieved  by  utilizing  the  more-complete- 
expansion  cycle.   If  the  engine  is  operated  frequently, 
the  amount  of  money  saved  in  fuel  is  significant  and 
could  pay  off  the  cost  of  the  engine  modification  in 
less  than  a  year.   Therefore,  this  idea  is  recommended 
for  further  development  and  testing. 


